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INTRODUCTION 


The  purpose  of  this  contract  was  to  study  beam  expansion  for  multimode 
fiber  optics  and  optical  signal  processing  applications.  It  was  proposed  to 
integrate  the  beam  expander  monolithically  with  a  DBR  laser  to  provide  an 
expanded  beam  about  75  pm  wide.  Simultaneous  studies  were  begun  to  develop 
the  beam  expander  to  demonstrate  the  concept  and  also  to  fabricate  DBR 
lasers.  It  was  the  intention  of  the  program  to  integrate  these  two  devices 
together  during  the  second  year  of  the  contract.  Within  the  first  year,  beam 
expansion  was  demonstrated  using  a  distributed  Bragg  deflector  (DBD)  in 
passive  material,  and  its  unique  polarization  properties  were  characterized. 
However,  unexpected  difficulties  were  encountered  in  the  fabrication  of  DBR 
lasers,  and  that  part  of  the  project  was  delayed  by  a  year  before  a  successful 
fabrication  schedule  was  worked  out.  The  details  of  the  fabrication  studies 
which  went  into  the  successful  fabrication  of  DBR  lasers  is  contained  in 
Aerospace  Report  No.  ATR-79(7739)-l,  submitted  by  The  Aerospace  Corporation  to 
University  of  Southern  California  on  April  10,  1979. 

During  the  second  year  DBR  lasers  were  successfully  fabricated  and  their 
properties  were  studied,  in  order  to  obtain  information  as  to  how  successfully 
they  would  integrate  with  the  beam  expander.  However,  the  contract  ran  out 
before  the  monolithic  integration  was  achieved. 

Because  of  the  delay  In  DBR  laser  fabrication  during  the  first  year  the 
program  was  broadened  to  include  studies  of  bistable  switches  for  fiber  optics 
and  optical  signal  processing  applications.  This  program  was  very  successful 
and  led  to  five  papers  and  four  talks,  listed  in  the  tables  at  the  end  of  the 
Introduction.  This  work  was  described  in  the  RADC  interim  report  (RADC-TR-79- 
112),  and  will  only  be  summarized  below.  During  the  second  year,  the 
sponsorship  of  this  work  was  picked  up  by  AFOSR,  and  further  research  in  this 
area  is  described  in  the  interim  scientific  report  submitted  to  AFOSR. 

The  realization  of  the  polarization  properties  of  the  distributed  Bragg 
deflector  beam  expander,  which  make  it  the  integrated  optics  analog  of  the 
Glan-Thomp8on  prism,  suggested  the  possibility  of  using  the  DBD  along  with 
polarization  modulation  to  made  a  nultlmode  fiber  optics  switch.  Work  was 
therefore  begun  on  developing  LiNb03  waveguides  with  good  optical  properties 
to  be  used  with  the  beam  expander  to  make  a  multimode  fiber  optics  switch. 


Each  of  the  efforts  just  described  will  be  summarized  in  the  following 
sections  of  the  introduction  and  then  will  be  described  in  separate  chapters 
in  the  body  of  this  report. 

1.  Waveguide  Beam  Expansion 

In  accordance  with  the  proposal,  studies  of  waveguide  beam  expansion 
were  undertaken  using  a  distributed  Bragg  deflector  on  VPE  GaAs  wave¬ 

guides  along  with  excitation  by  a  1.15  pm  HeNe  laser.  The  experiments 
resulted  in  a  deflected  beam  with  a  width  of  1300  pm,  for  an  expansion 
factor  of  26.  This  expansion  was  too  large  for  the  applications  intended 
in  the  contract,  however.  The  reason  for  the  large  expansion  was  the  fact 
that  the  grating  had  a  long  coupling  length.  SEM  photographs  of  the 
grating  showed  it  to  be  of  excellent  quality  but  since  the  VPE  waveguides 
had  a  large  effective  waveguide  thickness,  the  coupling  coefficient  was 
weak.  The  result  of  this  study  was  the  prediction  that  the  distributed 
Bragg  deflector  should  provide  adequate  beam  expansion  when  double 

heterostructure  material  is  used. 

The  experiments  demonstrated  that  the  TE  mode  was  not  deflected  when 
the  grating  spacing  and  angle  were  chosen  to  produce  a  90°  deflection. 

Only  the  TM  mode  was  deflected,  while  the  TE  mode  was  transmitted  through 
the  DBD.  Measurements  indicated  an  extinction  ratio  of  343  between  the 
transmitted  components.  TMs  suggested  the  use  of  the  DBD  as  a 

polarization  analyzer,  or  as  the  integrated  optics  analog  of  the  Glan- 

Thompson  prism.  However,  it  also  indicated  that  the  monolithically 
integrated  beam  expander  would  have  to  operate  at  45°,  rather  than  at  90°, 
since  the  DBR  laser  operates  with  the  TE  polarization. 

At  45°  beam  expander  was  fabricated  in  VPE  GaAs,  although  again  the 
expanded  beam  was  too  large  to  be  useful.  The  decision  to  use  the  45° 
deflection  angle  led  to  a  redesign  of  the  beam  expander,  which  could  now 
use  a  first  order  grating.  Measurements  were  made  of  the  coupling 

coefficient  of  the  gratings  used  in  the  DBR  lasers  and  led  to  a 

prediction  that  expansion  factors  of  5,  or  deflected  beams  the  order  of 
75  pm  should  be  easily  achieved.  The  result  of  this  study  is  the 
assurance  that,  with  the  redesign  to  45°  beam  deflection,  expansion 


factors  as  predicted  in  the  proposal  can  be  achieved  in  double 
heterostructure  material  which  we  have  available  at  the  present  time. 

DBR  Lasers  for  Fiber  Optics  and  Optical  Signal  Processing  Sources 

DBR  lasers  were  successfully  fabricated  during  the  second  year,  and  a 
number  of  studies  were  undertaken  to  evaluate  their  suitability  as  sources 
for  multimode  fiber  optics  and  optical  signal  processing.  We  observed  the 
single  mode  operation  characteristic  of  the  DBR  lasers,  but  also  observed 
that  if  the  wavequide  layer  was  too  thick,  two  wavelengths  appeared, 
separated  by  40  A.  This  is  due  to  the  fact  that  the  waveguide  supports 
two  modes  which  have  different  propagation  constants  and  therefore  satisfy 
the  Bragg  condition  at  different  wavelengths. 

We  made  measurements  of  the  f ront-to-rear  power  ratio  emitted  from 
the  lasers  in  order  to  determine  the  coupling  coefficient  of  the  gratings 
causing  retroref lection.  The  effective  1/e  coupling  lengths  of  the 
gratings  were  192  pm  and  106  pm,  numbers  which  were  in  good  agreement 
with  estimations  made  from  the  tooth  shape  and  heights  and  the  waveguide 
effective  thickness.  These  studies  confirmed  that  the  double 
heterostructure  material  could  be  used  to  make  excellent  beam  expanders. 

Finally,  we  measured  the  in-plane  beam  divergence  of  several  DBR 
lasers  and  found  one  whose  behavior  was  diffraction  limited,  right  above 
threshold,  and  would  be  suitable  for  a  beam  expander.  The  other  device, 
however,  displayed  multimode  characteristics  which  would  not  be  suitable 
for  controlled  beam  expansion.  This  suggests  that  it  will  be  necessary  to 
control  the  spatial  modes  of  the  DBR  laser  in  the  same  fashion  as  has  been 
recently  developed  for  Fabry-Perot  lasers,  using  narrow  stripes  and 
tapered  junctions,  if  optimum  performance  is  desired. 

Bistable  Optical  Switches  for  Multimode  Fiber  Optics  and  Optical  Signal 
Processing  Applications 

When  the  output  of  any  modulator  or  switch  is  detected  and  used  to 
provide  a  voltage  across  that  same  modulator  or  switch,  that  device  is 
bistable.  This  means  that  there  is  a  regime  of  operation  in  which  there 
are  two  stable  final  states,  depending  on  the  history  of  operation  of  the 
device  (whether  the  steady  state  was  reached  with  an  incident  signal  below 


or  above  its  steady-state  value).  We  demonstrated  this  concept  in  a  bulk 
modulator  and  suggested  that  an  integrated  optics  version  with  a  multimode 
bistable  optical  switch  could  be  used  for  multimode  fiber  optics  and 
optical  signal  processing  applications.  For  this  reason  we  embarked  on  a 
study  of  these  devices  in  the  first  year  of  this  contract. 

We  demonstrated  an  integrated  optics  version  of  the  bistable  switch 
using  a  multimode  out-diffused  LiNbO^  modulator  which  was  only  2  mm 
long.  Bistability  was  achieved,  and  the  details  are  described  in  the 
pubications  listed  in  Table  I  at  the  end  of  this  introduction. 

In  addition,  we  made  studies  of  the  transient  response  of  these 
bistable  optical  switches.  In  particular,  we  demonstrated  the  phenomenon 
of  "critical  slowing  down".  This  is  a  delay  in  the  switching  of  the 
device  which  occurs  if  the  incident  signal  is  near  its  critical  value  for 
switching.  This  phenomenon  means  that  the  switch  signal  must  be  bigger 
than  its  minimum  value  if  long  delays  are  to  be  avoided.  It  alsc  means 
that  estimates  of  the  required  switch  energy  based  on  steady-state  values 
of  incident  signal  and  time  constants  will  be  wrong.  Detailed 

understanding  of  this  phenomenon  and  presentation  of  the  experimental 
results  are  In  the  publications  listed  at  the  end  of  the  Introduction. 

In  the  past  year,  when  support  of  the  bistable  switch  as  picked  up  by 
AFOSR,  we  have  made  progress  in  further  understanding  the  transient 
behavior  of  these  devices.  In  particular,  we  have  studied  the  effects  of 
overshoots  which  occur  when  the  device  has  two  time  constants.  In 
addition  we  have  studied  the  distributed  feedback  non-linear  optical 
device,  which  can  also  be  bistable.  These  studies  are  described  in  the 
interim  report  to  AFOSR. 

4.  Studies  of  LiNbO^  Waveguide  Loss 

Motivated  by  the  possibility  of  Integrated  optics  multimode  switches 
in  LiNbOj  which  could  be  butt-coupled  to  multimode  fiber  optics,  and  by 
the  interest  in  LiNbOj  bistable  devices  for  optical  signal  processing, 
studies  were  undertaken  to  improve  the  characteristics  of  LiNbO^ 
waveguides.  In  Chapter  3,  studies  of  the  loss  mechanisms  in  LiNbO^ 
waveguides  are  described.  The  loss  contributions  are  separated  into 
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absorption,  scattering  from  the  surface  of  the  waveguide  and  scattering 
from  within  the  waveguide. 

Absorption  losses  were  measured  calorimetrically  and  shown  to  be 
negligible.  A  theoretical  study  was  undertaken  which  demonstrated  that 
waveguides  buried  below  the  surface  of  the  wafer  should  have  considerably 
less  scattering  loss  than  those  which  are  at  the  surface. 

Experimental  studies  indicated  that  scattering  from  within  the  bulk 
of  in-diffused  waveguides  was  much  smaller  than  the  scattering  from 
surface  roughness.  Finally,  insertion  loss  studies  were  made  to  compare 
the  relative  importance  of  length-dependent  and  input  coupling  losses. 
The  results  show  that  the  length-dependent  losses  in  the  waveguides  we 
fabricated  varied  from  0.6  dB/cm  in  waveguides  20  pm  wide  to  1.1  dB/cm 
in  waveguides  10  pm  wide  to  1.6  dB/cm  in  waveguides  3  pm  wide.  The 
input  coupling  loss  was  found  to  be  between  0.6  and  2  dB. 

Several  experiments  were  undertaken  in  an  attempt  to  reduce  the 
insertion  loss  of  LiNbOj  butt-coupled  multimode  waveguides.  These  include 
polishing  the  surface  and  applying  a  high  index  cladding  layer  to  the 
surface  of  the  waveguide.  The  experiments  did  not  reveal  any  sizeable 
improvement  in  waveguide  performance  as  a  result  of  these  experiments. 

5.  Studies  of  Multimode  Waveguide  Fabrication  in  c-propagating  LiNbO^ 

The  use  of  the  DBD  as  an  integrated  optics  Gian  Prism  suggested  a 
multimode  fiber  optics  switch  based  on  electo-optic  polarization 
modulation.  This  switch  can  work  in  LiNbOj  for  c-propagating  light,  when 
the  electric  field  is  parallel  to  the  plane  of  the  waveguide  (as  will 
occur  for  electrodes  on  either  side  of  the  channel  waveguide).  The 
geometry  of  this  switch  will  be  described  in  Chapter  4. 

With  the  switch  as  motivation  we  began  to  study  the  fabrication  of 
multimode  waveguides  in  c-propagating  LiNbO^.  We  investigated  the  use  of 
an  electric  field  during  the  diffusion  to  assist  in  obtaining  deep 
diffusions.  The  results  were  somewhat  inconclusive,  since  the  samples  had 
a  tendency  to  crack.  Experiments  were  still  underway  at  the  end  of  the 
contract.  The  data  in  Chapter  4  represents  a  progress  report  to  date, 
rather  than  a  definitive  estimation  as  to  whether  or  not  electro-thermal 
diffusion  will  be  useful  in  the  fabrication  of  multimode  waveguides. 
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CHAPTER  1.  WAVEOUTPF  FAV  EXPANSION 

This  program  was  begun  to  study  waveguide  beam  expansion  for  use  in  fiber 
optics  and  signal  processing  applications.  Beam  expansion  can  be  achieved  by 
using  a  distributed  Bragg  deflector  (DBD).  This  device  consists  of  a 
corrugated  waveguide,  whose  grating  is  slanted  at  an  angle  6  with  respect  to 
the  incident  beam.  The  geometry  of  the  DBD  is  shown  in  Fig.  1.  The  incident 
beam  enters  from  the  left  into  a  confined  slab  guide  with  a  corrugated 
surface.  Bragg  reflection  occurs  from  the  grating,  reflecting  the  beam  at  an 
angle  equal  to  the  incidence  angle.  However,  the  reflections  from  the 
corrugations  add  up  coherently  only  if  the  Bragg  condition  is  satisfied;  that 
is. 


*  =  2® _ 

2n  sin  9* 


(1-D 


In  addition  to  beam  expansion,  applications  of  the  DBD  include  beam 
splitters  (power  division),  deflectors  (as,  for  example,  to  create  ring 
geometries),  polarizers  and  analyzers  (with  90°  deflection)  and 
multiplexers/demultiplexers  (with  chirped  gratings). 

The  theoretical  analysis  of  the  device  was  completed  during  the  initial 
phases  of  the  contract  and  is  presented  in  a  paper,  entitled,  "Distributed 
Bragg  Deflector:  A  Multifunctional  Integrated  Optical  Device",  H.M.  Stoll, 

Applied  Optics,  _1_7_  2562  (1978).  For  this  reason,  only  the  results  of  that 
theoretical  analysis  are  presented  here. 

It  was  found  that  the  incident  light  beam  amplitude  is  reduced  while 
traveling  through  the  DBD  region  by  an  exponential  factor  with  distance, 
exp  (  -  yz),  due  to  deflection.  The  results  of  the  calculation  are: 


and 


yTM 


where 


'TE 


2 

cos  9r 


n2 

log  (— )  (co89b-  1]  -  cos9g 


2 


(1-2) 


(1-3) 
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f 


r  / 


/  (  2  2^3 
(_n2  -  n^jd 

3  ( 2l  -  1)  n,  X  t3 
z  o  e 


(1-4) 


The  quantities  are  as  defined  in  Fig.  1-1,  with  Xq  the  wavelength  of 

incident  light,  l  the  Bragg  scattering  order  and  te,  the  effective  guide 

thickness.  If  the  mode  is  near  cutoff,  t  =  t  +  1/6,  where  6  is  the 

e 

exponential  decay  profile  of  the  mode  in  the  substrate. 


Consider  the  deflection  of  a  TM  polarized  beam  at  90°  which  is  the 
configuration  of  greatest  interest  for  beam  expander  applications: 


1 


log  n2 


(1-5) 


We  have  made  use  of  the  fact  that  nj  =  air  “  1. 

It  will  be  convenient  to  relate  these  results  to  those  of  retroref lection 
in  a  corrugated  waveguide.  Coupled  mode  theory  predicts  that  a  corrugated 
waveguide  will  couple  forward  and  backward  waveguide  modes  with  a  coupling 
coefficient  K.  Light  originally  travelling  forward  transfers  to  the  backward 
wave  exponentially  with  length,  with  a  coefficient  given  by  K.  Thus  the 
effectiveness  of  a  grating  can  be  described  by  describing  its  coupling 
coefficient.  For  square  corrugations  such  as  those  used  to  derive  Eq.  1-3, 
the  coupling  coefficient  is 


K 


o  2  »  2  ,  : 

a  4n  fd  ! 

3  (2)1  -  1)  X  n  It  > 

g  e 


(1-6) 


A  related  parameter  is  the  characteristic  distance  at  which  the  power  in  the 
initial  mode  decreases  by  a  factor  of  1/e.  This  characteristic  length  is 
given  by  •  (2K)  *. 

We  then  write  the  DBD  exponential  coefficient  as 


4ng  108  ns 

*  (n2  -1) 

g 


(1-7) 


Since  the  deflected  beam  has  an  amplitude  which  decreases  with  distance 
along  the  DBD  according  to  an  exponential  with  a  coefficient  y,  the 
Intensity  width  of  the  expanded  beam  will  be  approximately  l/2y. 
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Figure  1-1  Distribution  Bragg  deflector:  (a)  top  view;  (b)  side  cross 
sectional  view;  and  (c)  end  cross-sectional  view.  Small  and  large 
(solid)  arrows  in  (a)  represent  incident  and  deflected  beams  respec¬ 
tively.  Dashed  sinusoid  in  (c)  represents  equivalent  core  field  dis¬ 
tribution  in  the  y  direction. 


(substrate) 


SLAB' 
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Defining  the  expansion  factor  as  the  width  of  the  deflected  beam  divided 
by  the  incident  beam  width,  we  write  the  expansion  factor  as 


F 


1 


1 

2YiW 


2 


-  y  -  ^ 

4n  log  n 
8  g 


2 


(1-8) 


where  we  have  defined  the  grating  by  its  characteristic  length. 

It  can  be  seen  that  it  is  easy  to  obtain  large  expansion  factors,  since 
the  characteristic  length  of  gratings  can  be  typically  larger  than  W. 

For  the  case  of  GaAs  waveguides,  ng  =  3.6,  and  f  =  0.56.  The  intention 
of  the  program  was  to  obtain  expansion  factors  of  roughly  five;  to  expand  the 
20  yin  emission  from  a  GaAs  stripe  geometry  laser  to  a  100  pm  beam,  which 
could  be  used  for  signal  processing  or  which  could  be  used  with  an  output 
grating  coupler  to  make  a  beam  which  could  be  butt-coupled  to  multimode 
fibers,  discussed  in  Chapter  2.  Expansion  of  a  factor  of  5  can  be  achieved  if 
L  /W  =  2.8.  For  an  incident  beam  of  20  pm,  this  requires  L  =  56  pm,  or  K 

o  o 

=*  .0089.  A  grating  with  this  coupling  coefficient  can  be  achieved  if  the 
waveguide  is  1  pm  thick  and  the  tooth  depth  Is  0.1  pm  (Eq.  1-6). 

Reflection  with  no  expansion  will  occur  when 

(n  2  -  1) 

W  =  sfl  - 1 -  L  =  \fl  f  L  .  (1-9) 

4  „  2  log  n  8  8 

g  6  8 

For  GaAs,  W  *  0.79  Lg  will  ensure  reflection  with  no  expansion,  at  90° 
deflection. 

The  research  program  was  designed  first  to  demonstrate  beam  expansion  in 
high  optical  quality  VPE  GaAs  waveguide  material  and  then  to  apply  the  concept 
to  LPE  GaAs  double  heterostructures,  so  that  the  laser  and  beam  expander  could 
be  combined  into  the  same  substrate. 

Demonstration  of  beam  expansion  in  VPE  GaAs  epilayers  was  successful.  A 

perspective  view  of  the  Bragg-effect  polarizer/analyzer  which  we  fabricated  is 

shown  in  Fig.  1-2.  The  asymmetric  slab  waveguide  structure  consisted  of  a 

✓  16  —3 

1  pm  thick,  n-type  (~10  cm  )  GaAs  layer  epitaxially  grown  over  a 
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Figure  1-2  Perspective  view  of  the  Bragg-effect  polarizer /analyzer 
which  was  tested.  The  electric  field  vector  of  the  input  beam  is 
oriented  at  an  angle  6  to  the  waveguide  surface  normal. 


heavily  doped  n-type  ("  5  x  10  cm  )  GaAs  substrate;  it  supported  only  the 

lowest  order  one-dimenional,  T.E.  -  and  T.M.  -  polarized  modes.  Ninety  degree 

Bragg  scattering  of  the  incident  1.15  pm  radiation  was  accomplished  via  a 

second  order  periodic  surface  corrugation  (A  **  0.478  ym)  which  was  ion  milled 

into  the  waveguide  surface  and  overcoated  with  a  0. 1  ym  thick  layer  of 
2 

aluminum  (n^  «  -69)  in  order  to  increase  its  deflection  efficiency.  As 
revealed  by  SEM  examination  of  a  cross  section  through  the  device,  the 
corrugations  bore  a  nearly  triangular  tooth  profile.  A  photograph  of  a 
typical  grating  is  shown  in  Fig.  1-3.  The  total  corrugation  length,  L,  and 
input  beam  width,  w,  were  3.2  mm  and  50  ym,  respectively. 

Figure  1-4  shows  a  semi-logarithmic  plot  of  the  near-field,  deflected 
beam  intensity  profile  of  one  of  our  90°  Bragg  deflectors.  The  data  shows 
that  the  deflected  beam  decreases  of  a  factor  of  1/e  in  1.3  mm  length.  With 
in  input  beam  of  50  ym,  this  corresponds  to  an  expansion  factor  of  38,  mich 
larger  than  required  for  fiber  optics  applications.  From  this  data  it  can  be 
inferred  that  the  effective  coupling  length  of  the  grating  was  L  =  385  ym, 
with  a  coupling  coefficient  of  K  =  0.0013.  This  coupling  is  too  weak  for 
fiber  optics  applications.  The  reason  for  the  weak  effect  is  use  of  VPE 
GaAs.  The  index  change  due  to  the  reduction  of  the  number  of  free  carriers  in 
the  epilayer  is  very  small.  As  a  result,  the  guided  mode  extends  into  the 
substrate  and  the  effective  guide  width  is  larger  than  the  1  ym  physical 
thickness  of  the  epilayer.  For  the  physical  parameters  used  in  the 
experiment,  the  effective  guide  thickness  was  several  microns,  causing  the 
coupling  coefficient  to  be  small. 

Because  of  the  weakness  of  guiding  in  the  VPE  material,  the  effective 
waveguide  thickness  was  large  and  the  grating  coupling  coefficient  was 
small.  As  a  result,  the  expanded  beam  was  too  large  to  be  useful  for  more 
than  demonstration  of  the  principles.  However,  it  was  demonstrated  that  even 
a  grating  with  a  low  coupling  coefficient  can  be  used  successfully  as  a 
polarizer-analyzer.  A  study  of  this  effect  was  undertaken  and  the  results  of 
this  study  have  been  submitted  for  publication  (H.M.  Stoll,  W.E.  Soady  "Bragg- 
effect  Polarizer/Analyzers  and  Lens-less  Waveguide  Beam  Expansion")  and  are 
summarized  in  this  report. 

Operation  of  the  DBD  as  a  polarizer/analyzer  can  be  seen  by  observing 
that  YfE  *  0  for  a  deflection  angle  of  90°.  Polarizer/analyzer  performance 
parameters  of  interest  include  the  extinction  ratio,  p,  which  we  define  as 
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-3  SEM  photograph  of  cross-section  of  ion-milled  grating 


being  equal  to  the  ratio  of  T.M.  -  to  T.E.  -mode  power  deflected  by  the 

device,  or 

yTM  6TE  1  “  exp  (-  2  6^  L)  , ,  ,  ns 

^TE  *  ^TM  1  -  e*P  2  «TE  L>  ’ 

and  the  T.M.  deflection  efficiency,  n,  which  we  define  as  being  equal  to  the 

ratio  of  deflected  to  incident  T.M.  -mode  power,  or 


■y 

n  -  ■—  [1  "  exp(-  2  5^  L)]  .  (1-11) 

TM 

In  Eqs.(l-lO)  and  (1-11),  6TE  and  6^,^  are  the  total  (including 
background)  T.E.  -  and  T.M.  -mode  attenuation  coefficients,  respectively,  and 
L  is  the  device  length. 

Experiments  were  performed  on  the  polarization  dependence  of  beam 
deflection  using  the  VPE  GaAs  waveguide  beam  expander  shown  in  Fig.  1-2.  By 
measuring  the  near  field  intensity  profile  of  the  beam  which  was  Bragg 
deflected  by  the  polarizer/analyzer,  6^  was  determined  to  be  4.19  cm 
Knowledge  of  6^  together  with  measurements  of  the  undeflected  throughput 
power, 

P  (0)  “  P  in  [c°s2e  exp('2  5TM  L)  (1-12) 

out  „ 

+  sin  0  exp(-2  6TE  L)] 

for  0  ■  0  and  ir/2,  were  then  used  to  yield  6TE"  3.11  cm* ,  ■  0.154  cm  * , 

and  n  ”  0.034.  Finally,  direct  measurement  yielded  an  extinction  ratio  of 

-4  -1 

343  from  which  it  was  inferred  that  6TE  ■  3.6  x  10  cm  .  Figure  1-5, 

which  plots  P  fc(9),  PTUO(0)  (the  total  T.M.  power  deflected),  and  the 
out  imo 

experimental  data,  summarizes  the  results. 

The  passive  beam  expander  in  VPE  GaAs  demonstrated  that  beam  expansion 
was  Indeed  possible  using  the  Distributed  Bragg  Deflector.  However,  it  was 
also  demonstrated  that  the  efficiency  of  deflection  goes  to  zero  for  the  TE 
mode  at  a  deflection  angle  of  90°.  Before  we  knew  whether  or  not  the  90° 
angle  of  deflection  was  suitable  for  use  in  a  monolithic  circuit  together  with 
the  laser  source,  we  needed  to  know  what  the  polarization  would  be  of  the  DBR 
lasers,  whose  characteristics  will  be  described  in  the  next  chapter.  It  was 
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found  that  these  lasers  operated  with  a  TE  polarization.  This  means  that  the 
angle  of  deflection  must  be  some  angle  other  than  90°.  The  most  convenient 
angle  of  deflection  is  45°. 

Figure  1-6  shows  a  45°  DBD  which  was  fabricated  in  VPE  GaAs.  The  beam 
expansion  is  shown  in  the  lower  photograph,  as  a  streak  to  the  left  of  the 
beam.  The  difficulty  in  observing  this  device  is  the  necessity  of  polishing 
the  end  face  to  obtain  the  output  at  45°.  The  advantage  is  that  the  grating 
can  be  used  in  fundamental  order.  The  deflection  of  the  DBD  in  the  VPE  wave¬ 
guide  was  very  weak  because  the  waveguide  was  so  thick.  Experiments  on  45° 
deflection  in  thinner  LPE  GaAs  waveguides  awaited  the  development  of  DBR 
lasers.  Because  the  fabrication  of  the  DBR  lasers  took  more  effort  than 
anticipated,  there  was  not  time  before  the  end  of  the  contract  to  test  the  45° 
beam  expander  in  double  heterostructure  material. 

It  is  possible  to  use  the  results  of  the  DBR  lasers  and  the  coupling 
coefficients  of  their  gratings  to  predict  the  behavior  for  the  45°  beam 
expander,  however. 


The  45°  beam  expander  used  in  conjunction  with  a  GaAs  DBR  laser  will  be 
used  with  the  grating  in  fundamental  order  and  a  spacing  of  0.338  pm.  To 

determine  the  expansion  factor  from  the  K  value  of  the  grating,  we  use  Eq.  (1- 
2)  and  obtain  y  =  (2K/tt)2  W  cos2QD  =  W  cos20d/(ii  L  )2.  We  wish  to 

1  t  D  D  g 

determine  whether  or  not  teeth  of  the  shape  used  for  the  DBR  laser  can  be 

satisfactory  to  produce  beam  expansion  of  the  requisite  factor  of  five.  First 
we  need  to  define  an  expansion  factor,  which  includes  the  fact  that  the 

deflection  is  at  an  angle  other  than  at  90°.  The  width  of  the  expanded  beam 

is  sin  9p/(2  Y-pg)'  This  gives  an  expansion  factor  of 


F 


8 


2 


sin  0g 

2 ~ 
COS  0fi 


(1-13) 


The  expansion  factor  increases  rapidly  as  9>  90°,  but  at  this  angle  the  TE 
reflectivity  *0,  so  there  is  no  inconsistency.  From  Eq.  1-6,  it  can  be  seen 
that  the  effective  length  of  a  grating  used  in  first  order  differs  from  the 
effective  length  of  a  grating  used  In  third  order  by  a  factor  of  3. 


The  DBR  lasers  discussed  in  Chapter  2  had  an  average  K  *  0.0035  and 

effective  grating  length  L  »  140  ym.  In  first  order  these  gratings  would 

have  L  *  47ym.  From  Eq .  1-13  it  can  be  seen  that  for  W~15  ym,  F  is  too 
8  6 
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e  1.6.  Geometry  for  45°  beam  expander  along  with  photographs  of  an  IR  television  monito 
field  image  of  guided  beam.  The  top  photograph  represents  the  unexpanded  beam  50  pm  wid 
e  proper  angular  alignment.  The  bottom  photograph  demonstrates  beam  expansion  with  the 
k  to  the  left. 


large  to  produce  a  beam  which  will  match  into  fibers.  Appropriate  beam 
expansion  will  require  larger  coupling  coefficients  K  and  smaller  L  . 

O 

Progress  was  made  during  the  second  year  in  developing  heterostructure 
material  with  thinner  waveguides,  so  that  K  could  be  made  smaller.  Wave¬ 
guides  of  half  the  thickness  were  fabricated,  but  the  contract  expired  before 
DBR  lasers  could  be  fabricated  or  beam  expansion  measured. 


CHAPTER  2  DBR  LASERS  FOR  FIBER  OPTICS  AND  OPTICAL  SIGNAL  PROCESSING  SOURCES 


It  was  proposed  to  fabricate  DBR  lasers  as  part  of  this  program  and  to 
use  these  as  sources  together  with  the  grating  beam  expander.  This  integrated 
laser  and  beam  expander  has  applications  both  as  a  source  which  can  be  butt- 
coupled  to  multimode  fiber  optics  and  also  for  optical  signal  processing.  The 
geometry  of  the  integrated  DBR  laser  and  45°  DBD  beam  expander  is  shown  in 
Fig.  2-1.  The  beam,  which  is  expanded  to  approximately  100  pm,  can  be 

deflected  acousto-optically  to  make  a  spectrum  analyzer,  or  tc  make  a  beam 
whose  output  is  one-dimensiona lly  steerable,  for  use  in  optic'll  signal 
processing.  The  application  for  multimode  fibers  suggested  expansion  to 
60  pm  and  including  a  grating  output  coupler  so  that  a  fiber  could  be  butt- 
coupled  to  the  device.  This  geometry  is  shown  in  Fig.  2-2. 

The  program  succeeded  in  fabricating  DBR  lasers,  but  did  not  produce  the 
integration  of  the  beam  expander  and  the  DBR  laser.  This  is  because  it  was 
discovered  that  it  took  longer  to  develop  techniques  for  fabrication  of  DBR 
lasers  than  originally  anticipated,  and  the  first  year  was  spent  in  developing 
a  processing  schedule  for  DBR  lasers.  This  processing  schedule  is  included  as 
an  appendix  to  this  chapter.  With  the  development  of  the  processing  schedule, 
it  was  possible  to  fabricate  DBR  lasers  during  the  second  year.  This  chapter 
will  summarize  the  measurements  made  on  those  lasers,  and  describe  how  these 
properties  affect  expected  performance  of  the  integrated  sources  for  fiber 
optic  and  optical  signal  processing. 

The  geometry  of  the  DBR  lasers  which  were  fabricated  is  shown  in  Fig.  2- 

3.  The  gratings  were  third  order,  and  the  laser  region  was  a  large  optical 

cavity.  The  design  studies  for  this  configuration  were  given  in  the  interim 

report.  This  design  was  not  optimum  because  the  waveguide  region  was  thicker 

than  desired,  but  was  the  best  material  available  at  the  time. 

Measurements  were  made  of  the  f ront-to-rear  power  ratio  of  DBR  lasers  of 
differing  lengths  to  determine  the  coupling  coefficient  of  the  grating  and  to 
check  for  its  reproducibility.  Two  devices  were  compared  in  detail;  the  data 
is  shown  in  Table  2-1. 


Figure  2-1  Distributed  Bragg  Deflector,  acting  to  expand  the  beam 
from  a  DBR  laser  into  a  wide  beam  within  the  plane  of  the  waveguide 
for  use  in  optical  signal  processing. 
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Figure  2-2  Distributed  Bragg  Deflector  plus  grating  output  coupler 
used  to  shape  the  beam  from  a  stripe  geometry  laser  to  the  mode 
pattern  of  a  multimode  fiber.  The  fiber  attaches  to  the  surface 
of  the  GaAs  double  heterostructure  wafer,  in  a  fashion  analogous 
to  the  Burrus  diode. 


.  34 ,  A  =  0.3  6  um 


Table  2-1 


Data  to  determine  coupling  coefficient  of  gratings  In  DBR  lasers 

Sample  Number  g9  g5 

Power  out  grating  (mW)  1  0.5 

Power  out  cleave  (mW)  30  10 

Power  ratio  (grating  to  cleave),  T  0.03  0.05 

Ln  T/4  -5.3  -4.4 

Grating  Length  1000  pm  464  pm 

K  (pm_1)  0.0026  0.0047 

Lg  =  ( 2K)-1  192  pm  106  pm 


According  to  the  theory,  a  retro-reflector  of  length  L  on  resonance,  has 
the  following  behavior  for  transmission  and  refection: 


Transmission: 


cosh  KL 


(2-1) 


Reflection: 


slnh2KL 

cosh^KL 


(2-2) 


When  KL  »  1: 


T  -  4e"2KL,  R  =  1  -  T 


(2-3) 


The  coupling  coefficient  was  determined  from  Eq.  (2-3)  using  the  fact 
that  the  f  ront-to-back  power  ratio  will  be  the  same  as  the  ratio  of  the 
transmission  of  the  grating  to  the  transmission  of  the  GaAs  cleave  on  the 
front  output  face  (85%). 

The  results  give  an  average  coupling  coefficient  of  K  =  0.0036  um  1 . 
We  now  compare  these  measurements  with  theoretical  expectations. 

The  theoretical  expression  for  the  coupling  coefficient  of  a  grating  was 
given  for  square  teeth  in  Chapter  1.  The  experimental  grating  looked  more 
like  triangular  teeth  than  square  teeth,  however.  The  comparison  between 
square  and  triangular  shapes  may  be  made  by  comparing  the  fourier  transforms 
for  the  third  order  coefficients  for  square  and  triangular  shapes.  Using  the 
calculations  shown  in  Fig.  2-4,  we  see  that  the  polarization  induced  by  a 
triangular  grating  will  be  2/3  that  induced  by  the  step  index.  Thus  we  expect 
the  grating  coupling  coefficient  to  be 


2  2 

4t  (n  -1)  r  a_ 
9mXn  l  t 


[  1  + 


2-rr  V  n  -  1 


3  OJaT 
n2  -  1 


(2-4) 


(n  ■  index  of  guide) 


where  we  have  taken  the  result  of  Yariv,  multiplied  by  2/3  and  used  teff 
instead  of  t  for  guides  which  are  not  infinitely  well  confining. 
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Figure  2-4  Calculation  of  third  harmonic  for  square  and  trianeular 
shape  grating  teeth.  8 
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From  SEM  photographs  of  the  grating,  we  estimate  a  =  A/3,  t  =  1.54  pm, 

A  =  0.36pm.  The  grating  was  used  in  third  order,  and  the  refractive  index 

of  GaAs  at  0.88  pm  is  3.59.  These  numbers  give  the  theoretically  expected 

coupling  coefficient  of  K.  «*  0.0044pm  In  order  to  calculate  this  number, 

it  was  required  to  know  the  effective  guide  width.  This  was  determined  by 

solving  the  eigenvalue  equation  for  the  waveguide,  to  determine  the 

exponential  profile  outside  of  the  guide,  and  calculate  t  =  t  +  1/y.  For 

e  2  2 

the  DBR  lasers,  the  V  parameter  for  the  waveguide  was  AekQ  d  =  146, 
resulting  in  two  modes,  given  by  the  following  eigenvalue  equation  : 


tan  u 


1 


The  solutions  are  kd  =  1.7,  3.4,  and  give  yd  **  12.  This  makes  the  effective 
thickness  *  1.62  pm.  This  is  the  value  used  in  Eq.  (2-4)  to  determine  the 
theoretical  coupling  coefficient. 

The  excellent  agreement  between  theory  and  experiment  allows  us  to 
predict  behavior  for  other  waveguides.  Although  the  contract  ran  out  before 
any  more  lasers  were  studied,  we  have  fabricated  a  new  set  of  material  with 
waveguide  thickness  of  0.6  pm.  These  samples  should  yield  devices  with 
coupling  coefficients  as  much  as  eight  times  larger. 

Other  properties  relevant  to  DBR  lasers  as  sources  for  fiber  optics  and 
signal  processing  were  investigated.  These  include  the  spectrum,  total  power, 
the  P-I  characteristics,  and  the  beam  divergences  in  both  planes,  including 
whether  or  not  the  light  was  diffraction  limited.  The  results  are  summarized 
here.  A  typical  spectrum  of  a  DBR  laser  is  shown  in  Fig.  2-5.  First,  in 
order  to  prove  that  the  laser  action  occurred  from  the  grating  and  not  from 
spurious  reflections  the  following  tests  were  applied.  First,  it  was  observed 
that  the  spectrum  was  a  single  mode.  Comparable  spectra  for  Fabry-Perot 
lasers  from  the  same  material  were  multimode,  as  shown  in  Fig.  2-6.  Secondly, 
a  DBR  laser  which  operated  at  room  temperature  was  cooled  to  liquid  nitrogen 
temperature,  and  laser  action  ceased,  since  the  reflectivity  of  the  grating 
was  no  longer  at  the  frequency  of  the  fluorescence.  This  is  in  contrast  to 
Fabry-Perot  lasers  whose  frequency  shifts  with  temperature. 
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Figure  2-5.  Spectrum  of  typical  DBR  laser. 


(nm) 


In  some  of  the  diodes  two  lines  occurred,  as  shown  in  Fig.  2-7.  This 

occurred  because  the  waveguide  was  sufficiently  thick,  that  not  only  the  lowest 

even  mode,  but  also  the  lowest  odd  mode  could  be  confined  in  the  passive 
region.  These  modes  are  shown  in  Fig.  2-8.  These  two  modes  occur  in  the 

following  way. 

The  laser  wavelength  is  determined  by  the  condition  that  the  corrugation 
spacing  is  in  resonance  with  the  longitudinal  propagation  vector:  i.e. 

k^  =  3  rr/A.  However,  the  different  modes  have  different  propagation 

constants  since  the  transverse  propagation  constant  is  different  for  each 

mode.  The  dispersion  relation  thus  becomes 
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(m  +  1 )  tt 
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(2-5) 


To  get  a  simple  expression  for  the  difference  in  frequencies  for  DBR 
laser  lines  operating  on  different  waveguide  modes,  assume  that  the  effective 
waveguide  thickness  is  independent  of  mode,  that  the  refractive  index  does  not 
depend  on  frequency,  and  that  the  difference  in  wavelength  between  the  two 
laser  modes  Is  much  less  than  the  wavelength.  Then  the  equation  for  laser 
wavelength,  given  by 


4n  2 
— 2_ 
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- 

eff 


(2-6) 


can  be  differentiated,  yielding, 
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1.6  pm,  the  difference  in  wavelength 


For  an  effective  wave  guide  thickness  of 
between  modes  should  be  75A.  In  devices  tested,  the  actual  difference 

between  modes  was  40A,  indicating  that  the  assumption  that  the  waveguide 
refractive  index  was  a  constant  value  was  only  rough  approximation.  More 
exact  results  can  be  obtained  by  solving  numerically  for  the  eigenmodes  of  the 
waveguide  and  numerically  solving  the  above  equation.  This  has  been  done 
graphically  in  Fig  2-9.  It  can  be  seen  that  the  laser  is  expected  to  operate 
at  wavelengths  in  good  agreement  with  the  experimental  results.  Finally,  it 
is  observed  that  the  longer  wavelength  is  the  lower  order  mode,  and  therefore 
the  one  with  the  lowest  threshold,  in  agreement  with  experiment. 
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The  linewidths  just  above  laser  threshold  were  typically  about  1A. 
Diode  G8,  sihown  in  Fig.  2-5,  Had  a  width  of  about  1.1  and  diode  G9,  shown 

in  Fig.  2-10,  had  a  width  of  1.6  A.  The  apparent  linewidth  broadened 
slightly  with  pumping,  as  shown  in  Fig.  2-10,  in  which  diode  g5  is  pumped 
three  times  threshold.  The  line  had  broadened  to  about  3  A,  but  is  was 
observed  on  the  oscilloscope  that  the  frequency  hanged  during  the  100  nsec  of 
the  input  pulse.  The  broadening  was  a  chirp  effect  due  to  heating  in  the 
junction.  Such  effects  have  been  seen  before  in  Fabry-Perot  diodes,  but  here 
are  much  smaller  because  of  the  smaller  temperature  sensitivity  of  DBR  lasers. 

Diode  g9  had  a  complicated  behavior,  since  the  two  frequencies 
corresponding  to  two  guided  modes  occurred  at  high  pumping  rates.  It  was 
observed  that  at  low  excitation  levels  the  shorter  wavelength  was  the 
strongest  peak  and  the  longer  wavelength  occurred  only  during  the  initial 
portion  of  the  pulse.  Evidently  heating  caused  the  resonance  condition  to 
change  and  near  the  end  of  the  pulse  the  shorter  wavelength  had  a  lower 
threshold  and  was  the  only  frequency  above  threshold.  This  data  demonstrates 
that  even  when  a  grating  is  used,  it  is  possible  to  get  several  wavelengths 
oscillating,  not  necessarily  simultaneously  in  time  within  the  laser.  Clearly 
DBR  diodes  must  be  optimized  to  have  well-behaved  results. 

The  complexity  of  this  device  can  be  seen  from  the  P-T  curve  (Fig.  2-11), 
which  are  filled  with  kinks.  The  DBR  laser  causes  only  one  or  two  frequencies 
to  oscillate,  but  does  not  necessarily  improve  P-I  characteristics  over  Fabry- 
Perot  lasers. 

Since  kinks  seem  to  be  related  to  the  filamentary  structure  of  lasers, 
and  since  we  were  interest  in  beam  quality,  in  order  to  perform  beam 
expansion,  we  made  a  careful  study  of  the  diffraction  qualities  of  the  light 
emitted  from  the  DBR  lasers.  Depending  on  the  device,  we  were  able  to  observe 
diffraction  limited  operation,  just  above  threshold.  This  is  demonstrated  by 
studying  the  divergence  of  sample  g9,  right  above  threshold.  Near  field 
patterns  were  observed  both  out  the  cleave  and  the  grating  face  (Fig.  2-12). 
In  addition,  the  in-plane  beam  divergence  was  measured  coming  out  of  the 
cleave.  This  was  done  by  monitoring  the  spot  size  as  a  function  of  distance 
from  the  cleave.  The  data  are  shown  in  Fig.  2-13. 
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Figure  2-10  Spectra  of  diodes. 
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Figure  2-12  Near -field  profiles  of  DBR  laser.  Diode  g9,  at  threshold 


BEAM  WIDTH,  fim 


Figure  2-13  In -plane  beam  width  as  a  function  of  distance  from  the 
exit  face  (as  measured  on  a  television  monitor).  Far  field  divergence 
angle  is  72  m  rad. 


The  data  are  analyzed  by  comparing  the  measured  divergence  with  that 
expected  for  a  Gaussian  beam  of  the  width  determined  from  the  near  field 
intensity  profiles  (Fig.  2-12). 

The  angular  divergence  of  a  Gaussian  beam  is  given  by 

9  =  X / 7i  a)  (2-8) 

o  o 

'  2 

where  9  is  the  half-angle  and  u  is  the  half  width  of  the  beam  at  1/e  in 
o  o 

intensity. 


Redefining  the  quantities  in  terms  of  full  width  half  maximum, 
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X _ 

0. 72iru) 


(2-9) 


The  filament  measured  in  the  near  field  had  an  apparent  width  of 
6.4  pm  (correction  has  not  been  made  for  resolution  of  the  microscope 
objective  used  in  this  measurement).  The  far  field  beam  divergence  out  the 
cleave  was  0.072  rad,  which  corresponds  to  the  diffraction  expected  from  a 
filament  of  width  5.4  pm.  This  is  in  good  agreement  with  the  measured 
filament  width,  and  supports  the  assumption  that  this  laser  operated  in  a 
single  spatial  mode. 

The  laser  output  from  the  grating  face  was  observed  to  have  a  near  field 
profile  which  was  a  single  filament  22  pm  wide.  This  profile  is  the  beam 
which  has  diffracted  from  the  5.4  pm  filament  in  the  laser  region  due  to  the 
use  of  a  1  mm  long  grating  as  shown  in  Fig.  2-14.  If  the  filament  is  assumed 
to  extend  the  length  of  the  active  region  of  the  laser  diode,  and  begins  to 
expand  when  it  enters  the  grating,  then  the  diffraction  should  occur  for  a 
distance  of  1  mm,  at  an  angle  given  by  the  far  field  divergence  angle  divided 
by  the  index  of  refraction  of  the  medium  (0.02  rad).  This  divergence  would 
cause  the  output  of  the  grating  to  have  a  width  of  20  pm.  This  compares 
favorably  with  the  measured  filament  width  of  22  pm. 

For  completeness,  out-of-plane  beam  divergence  was  measured  also.  It  was 
observed  that  the  frequency  corresponding  to  the  higher  order  waveguide  mode 
caused  a  double-humped  far  field  pattern  from  the  cleave.  The  divergence  was 
greater  than  0.48  rad  and  could  not  be  measured  since  it  was  limited  by  the 
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Figure  2-15  Near  field  profiles  of  DBR  laser.  Diode  g5. 
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aperture  of  the  lens  which  was  used.  The  other  frequency  was  a  lowest  order 
waveguide  mode  and  had  a  beam  divergence  of  0.3  rad.  This  proves  the 

assumption  discussed  earlier  that  the  two  frequencies  are  due  to  different 
modes  occurring  in  the  laser  region. 

Other  diodes  did  not  necessarily  show  in-plane  diffraction  limited 

behavior.  The  near  field  profiles  of  diode  g5  are  shown  in  Fig.  2-15. 

Although  the  apparent  near  field  profile  was  15  pm  wide,  the  beam  divergence 
was  0.16  rad,  corresponding  to  a  filament  size  of  2.4  pm.  In  other  words, 
the  diode  had  a  highly  filamentary  character  which  caused  a  great  deal  of  beam 
divergence.  This  fact  was  made  clear  by  looking  at  the  profile  out  of  the 
grating.  The  output  had  several  filaments  and  an  overall  width  of  almost 

100  pm.  Diodes  such  as  this  would  not  be  very  useful  for  signal  processing 
applications  because  of  large  beam  divergence. 

The  result  of  this  study  is  the  fact  that  although  DBR  lasers  control  the 
longitudinal  modes,  they  do  not  improve  the  spatial  modes  or  ensure 

diffraction  limited  operation.  It  is  necessary  to  control  the  spatial  modes 
in  order  to  obtain  diffraction  limited  guided  beams  which  can  be  expanded  and 
used  for  coupling  to  fiber  optics  and  for  signal  processing.  The  spatial 
control  must  be  obtained  by  the  same  techniques  used  in  producing  single  mode 

Fabry-Perot  lasers:  narrow  stripes,  tapered  junctions,  and  other  modes 

control  devices.  Experiments  to  add  these  devices  to  the  DBR  lasers  were 

beyond  the  scope  of  this  project. 


APPENDIX 

Process  Schedule  -  DBR  Laser  Fabrication 

I.  Cap  layer  Examination 

1.  Cleave  sample  from  main  specimen 

2.  Etch  sample  for  5  sec.  in  1:1:8  (A:B:DI  H20)  A-B  etch. 

3.  Photograph  all  four  corners  of  sample  to  determine  cap  thickness 

variations . 

II.  Cap  Layer  Removal 

4.  Pre-photoresist  application  sample  cleaning:  Swab  with  methanol; 
flow  trichloroethylene  over  surface;  blow  dry  with  dry  nitrogen. 

5.  Spin  coat  sample  with  (unthinned)  AZ1350J  photoresist  @  5000  RPM 

for  30  sec.  Allows  removal  of  up  to  54  microns  of  gallium 

arsenide. 

6.  Pre-exposure  age  coated  sample  at  room  temperature  for  30  min. 

7.  Expose  coated  sample  for  45  sec.  using  mask  consisting  of  0.5  mm 
wide  bars  separated  by  1.0  mm  clear  spaces. 

8.  Develop  photoresist  for  1  min.  using  1:1  developer  (Shipley 
MF312:DI  H20). 

9.  Post  exposure  bake  at  100°C  for  1  hr. 

10.  Ion  mill  to  remove  cap  layer:  0.931  mA  beam  current;  1  KV 

accelerating  potential;  4.6  x  10-^  atmos.  argon. 

III.  P-isolation  Layer  Removal 

11.  Observe  photoluminescence  of  valleys  in  order  to  be  assured  of 

proper  cap  layer  removal.  Check  mesas  for  inadvertent  photoresist 
removal. 

12.  Chemically  etch  p-isolation  layer  and  active  layer  aluminum 

concentrations. 

13.  Observe  photo luminesce nee  of  valleys  in  order  to  be  assured  of 

proper  p-isolation  layer  removal.  Check  mesas  for  photoresist 
integrity. 
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IV.  Active  Layer  Removal 

14.  Ion  mill  to  remove  active  layer:  0.931  mA  beam  current;  1  kV 

accelerating  potential;  4.6  x  10-^  atmos.  argon. 

15.  Observe  photoluminescence  of  valleys  in  order  to  be  assured  of 
proper  active  layer  removal. 

16.  Remove  photoresist  with  plasmaline  etcher  (using  00  for  10  min). 

V.  Stripe  Geometry  Mask  Application 

17.  Clean  sample  in  boiling  xylene:  isopropyl  alcohol  (1:1)  for  2-3 
minutes.  Blow  dry  with  dry  nitrogen. 

18.  Deposit  3000  A  thick  Si02  layer  over  entire  sample  surface. 

19.  Pre-photoresist  application  sample  cleaning:  Swab  with  methanol; 
flow  trichloroethylene  over  surface;  blow  dry  with  dry  nitrogen. 

20.  Spray  deposit  5000  A  layer  of  thinned  AZ1350B  (A7.1350B;  thinner  - 
1:1).  Spray  for  0.77  sec.  @  5  inches  from  sample. 

21.  Pre-exposure  age  coated  sample  at  room  temperature  for  30  min. 

22.  Expose  coated  sample  for  45  sec.  using  dark  room  exposure 
machine.  Mask  consists  of  a  four-row  matrix  of  25  pm  by  400  pm 
windows. 

23.  Develop  photoresist  for  1  min.  using  1:1  developer  (Shipley 
MF312:DI  H20). 

24.  Post  develop  bake  at  100°C  for  1  hr. 

25.  Etch  Si02  strip  windows  with  plasmaline  etcher  (using  CF^  for  20 
min. ). 

26.  Remove  photoresist  with  plasmaline  etcher  (using  02  for  20  min.). 

VI.  Laser  Metallization  and  Fabrication 

27.  Clean  sample  in  boiling  xylene:  acetone  (1:1)  for  3  min.  Rinse 
in  three  changes  of  hot  xylene.  Blow  dry  with  dry  nitrogen. 

28.  Mount  sample  on  plating  fixture  with  "Dotite"  paint.  Coat  all 
areas  except  those  to  be  plated,  or  those  areas  already  protected 
by  photoresist,  with  apiezon  solution  (apiezon  wax  in  xylene). 
Let  dry. 

29.  Immerse  sample  in  HC1:H20  (1:1)  for  3  min.  Do  not  rinse. 

Transfer  directly  to  plating  solution. 


t  / 
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30.  Platinum  plate  back,  contact:  20v  open  circuit  voltage,  100  mA. 

2 

Current  for  3  rain,  (for  ~0.25  cm  sample). 

31.  Repeat  step  28  for  front  contact. 

32.  Repeat  step  29. 

33.  Platinum  plate  front  contact:  20v  open  circuit  voltage,  50  mA 

current  for  2  min.  (insufficient  trials  have  been  made  to 

determine  optimum  current-time  parameters  for  pattern  area 

plating,  but  these  values  have  been  used.  For  Fabry-Perot  samples 

2 

of  approximately  0.25  cm  area  time  should  be  increased  to  3 
min. ) . 

34.  De-mount  sample.  Clean  with  hot  xylene  and  blow  dry  with  dry 
nitrogen. 

35.  Dice  samples 

36.  Mount  on  headers  and  bond. 

VII.  Test  Active/Passive  Fabry-Perot 

37.  Measure  area  of  p-side  contact  and  begin  pulse  testing  at  1000 

2 

Amps /cm  .  (Total  current  passing  through  diode,  I,  is  related  to 
monitor  voltage,  V,  by  I  =  V/0.27.). 

VIII.  DBR  Grating  Application 

38.  Pre-Photoresist  application  sample  cleaning:  Swab  with  methanol; 
flow  trichloroethylene  over  surface;  blow  dry  with  dry  nitrogen. 

39.  Spray  deposit  ~0.15  pm  thick  layer  of  thinned  AZ1350B  photoresist 
(AZ1350B:  thinner  1:1).  Spray  for  0.125  sec.  @  5”  from  sample. 

40.  Pre-exposure  bake  @  65°C  for  25  min. 

41.  Expose  with  interfering  0.4579  pm  Argon  laser  beams  for  35  sec  @ 
combined  beam  intensity  of  10mW/CM^  (0.370  VDC  on  volt  meter  as 
detected  by  solar  cell).  Sample  is  aligned  using  He-Ne  laser 
prior  to  exposure.  Null  standing  wave  pattern  observed  at  beam 
splitting  prism.  Both  laser  beams  must  be  of  equal  intensity. 

42.  Develop  photoresist  in  diluted,  Shipley  MF312  developer  (MF312:DI 
H20  -  1:1). 

43.  Measure  grating  periodicity  using  ellipsometer ,  littrow  angles  of 
left  and  right,  first-order  diffracted  beams  are  measured  to 
within  20  sec.  of  arc  and  averaged.  Grating  orientation  with 


4  7 

i 


\  . 


■  r 


t 


respect  to  crystal  cleave  Is  simultaneously  measured  using 
rotatable  stage  mounted  on  ellipsometer.  The  latter  measurement 
is  accurate  to  within  one  minute  of  arc. 

44.  If  either  the  grating  period  or  grating  orientation  is  wrong,  beam 
angle-of  incidence  must  be  adjusted  and  steps  38  through  43 
repeated. 

IX.  Test  DBR  Laser 

45.  Same  as  step  37. 
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CHAPTER  3.  STUDIES  OF  LiHb03  WAVEGUIDE  LOSS 

The  research  at  USC  during  the  first  year  centered  on  the 
development  of  the  multimode  bistable  optical  switch.  This 
work  was  described  fully  in  the  interim  report,  and  in 
several  publications  listed  in  the  introduction.  Continua¬ 
tions  of  this  research  were  funded  by  AFOSR  allowing  research 
on  this  contract  during  year  two,  to  focus  on  the  sources  of 
loss  in  LiNbO^  waveguides,  and  on  efforts  to  fabricate 
LiNbO^  multimode  polarization-independent  waveguides.  This 
work  was  motivated  by  the  fact  that  good  quality  waveguides 
and  also  multimode  polarization-independent  modulators  do  not 
exist,  and  because  there  are  no  low-loss  multimode  and 
polarization-insensitive  waveguides  in  LiNbO^.  The  first  part 
of  this  study  involved  investigations  of  losses  in  LiNbO^ 
waveguides.  The  second  part  of  the  study  involved  searching 
for  techniques  to  fabricate  LiNbO^  multimode  waveguides 
which  guide  both  TE  and  TM  modes  and  be  mode-matched  to 
multimode  fibers.  This  is  described  in  Chapter  4. 
a.  Sources  of  Loss  in  Waveguides 

In  waveguides  fabricated  by  a  diffusion  process  such  as 
LiNbO^,  mechanisms  for  loss  are  absorption,  scattering  within 
the  waveguide  and  scattering  from  the  surface  of  the  waveguide. 
Absorption  losses  may  come  from  either  the  bulk  substrate 
material  or  from  the  diffusant.  If  the  substrate  is  trans¬ 
parent  and  the  diffusant  is  properly  chosen,  absorption  loss 


should  be  negligible.  We  used  calorimetry  to  measure  the 
absorption  losses  in  LiNbO^  waveguides,  experiments  which  will 
be  discussed  in  Section  c  below.  We  found,  indeed,  that  with 
proper  preparation,  this  loss  mechanism  was  much  lower  than 
scattering  loss  in  LiNbO^  for  both  in-diffused  and  out-diffused 
waveguides . 

There  are  several  sources  of  scattering  loss.  The  first 
is  from  the  surface  of  the  waveguide.  Any  imperfections  in 
the  flatness  of  the  waveguides  causes  scattering  from  guided 
modes  into  radiation  modes.  Because  of  this  it  is  reasonable 
to  assume  that  the  loss  would  be  lower  for  waveguides  which 
are  buried  below  the  surface  of  the  waveguide.  In  order  to 
test  the  validity  of  this  idea,  we  calculated  the  power  loss 
due  to  surface  imperfections  for  waveguides  which  were 
diffused  at  the  surface  and  compared  them  to  those  which  were 
buried  below  the  surface.  We  compared  these  results  with 
the  losses  from  step-index  guides.  We  found,  as  expected, 
that  the  losses  were  lowest  for  guides  which  were  buried 
below  the  surface.  This  theoretical  study  is  described  in 
Section  b. 

Scattering  may  also  occur  from  nonuniformities  in  the 
refractive  index  caused  by  the  diffusant  or  present  in  the 
substrate  material.  These  nonuniformities  in  refractive 
index  may  be  due  to  aggregates  of  the  diffusant,  occlusions, 
crystallites,  or  variations  in  stoichiometry  or  composition. 
These  refractive  index  fluctuations  will  be  small,  as  long 
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as  the  bulk  substrate  material  is  of  optical  quality  and  as 
long  as  the  diffusant  enters  the  material  in  a  sufficiently 
dispersed  fashion.  In  order  to  investigate  the  amount  of 
scattering  from  the  nonuniformities  in  the  refractive  index 
and  compare  that  to  the  amount  of  scattering  from  the  surface 
of  the  waveguide,  we  made  comparative  studies  of  losses  and 
identified  the  mechanism  of  loss  in  LiNbO^  waveguides. 

These  experiments  are  described  in  section  d. 

As  part  of  our  experimental  studies  of  waveguide  losses, 
we  tried  several  experiments  to  reduce  the  insertion  loss 
for  LiNbOj  multimode  waveguides.  We  found  that  the  quality 
of  the  polish  on  the  end  faces  was  an  important  factor  in 
reducing  waveguide  losses.  Other  techniques  we  tried  to 
reduce  the  losses  were  not  so  satisfactory.  The  details  of 
these  experiments  are  described  in  Section  e. 
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b.  A  Comparison  of  Losses  in  Imperfect  Surface-Diffused  and 
Buried  Optical  Waveguides 

The  recent  advances  in  integrated  optics  technology  have 
put  a  great  deal  of  emphasis  on  LiNbO^  and  LiTaO^  devices 
whose  waveguides  are  fabricated  by  diffusion.”''  In  most  cases, 
these  waveguides  are  formed  near  the  surface  and  have  graded 
refractive  index  profiles  usually  represented  by  exponential 
functions.  Early  experimental  work  on  diffused  glass  wave¬ 
guides  indicated  that  buried  guides  had  considerably  lower 
loss  than  surface  guides  '  .  In  this  section  we  present  the 
results  of  a  theoretical  analysis  which  show  that  in  principle, 
the  loss  in  a  buried  guiding  layer  formed  by  diffusion  is  at 
least  one  order  of  magnitude  lower  than  that  of  a  diffused 
waveguide  in  which  the  guiding  layer  adjoins  the  surface. 

These  results  suggest,  therefore,  that  for  low  loss  applications, 

4 

LiNbO-j  and  LxTaO^  waveguides  should  be  buried  to  minimize 
scatter  losses. 

We  assume  single  mode  waveguides  and  consider  the  loss 
due  to  radiation  from  waveguide  surface  imperfections.  It  is 
the  intent  of  this  paper  to  consider  a  number  of  comparable 
waveguide  geometries  to  determine  which,  in  principle,  has  the 
lowest  loss.  For  this  purpose  it  is  sufficient  to  use  sinusoidal 
wall  perturbation. 

We  shall  compare  graded  index  guides,  formed  by  diffusion, 
with  step  index  guides,  such  as  formed  by  thin  film  deposition 
or  crystal  growth.  In  both  cases  we  compare  both  surface 
(asymmetric)  and  buried  (symmetric)  guides.  It  will  be 
shown  that  the  buried  guides  have  considerably  lower 


loss-  In  addition,  it  will  be  shown  that  graded  index  wave¬ 
guides  are  inherently  less  sensitive  to  short  wavelength 
surface  roughness  than  are  abrupt  uniform  guides. 

The  approximate  formula  for  the  power  loss  is  obtained 
by  assuming  that  the  fields  of  the  perfect  guide  are  only 
slightly  perturbed  by  wall  imperfections,  which  can  be 
shown  to  yield  a  coupled  system  of  differential  equations 
for  the  coefficients  of  the  perturbed  fields'*  ^ .  Discrete  solu¬ 
tions  are  then  obtained  with  the  help  of  orthogonality  relations, 
together  with  initial  and  final  conditions.  The  general  solution 


for  scattering  loss  in  a  single  mode  guide  is  given  by: 


-  kQ  {/*  dp/*  dz  r  E(p)  n(x,z) 

4«  y  P  o  J  J 


E  dx} 


o 


where  ECp')  is  the  radiation  field,  E  is  the  guided 

mode  field,  and  n  (x,z)  is  the  refractive  index  discontinuity 


function  defined  as: 


n(x-,z)  =  n  (x,z) 


perfect 


n  (x, z)  imperfect  (2) 


It  can  be  shown  that  the  resulting  calculation  performs  the 
fourier  transform  of  wall  imperfections.  Thus  the  physical 
processes  and  comparisons  of  different  geometries  of  waveguide 
profile  can  be  determined  by  looking  at  sinusoidal  perturba¬ 
tions 
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First  we  discuss  symmetric  buried  layers  whose  index 


profile  is  represented  by  the  following  distribution: 

n(x)  =  n  +  in  cosh  (x  -  x  )/d)  (3) 

s  o 

where  ng  is  the  substrate  index  and  An  is  the  peak  index 
discontinuity.  Solutions  of  the  guided  TE  modes  for  this 
profile  have  been  obtained  previously  in  the  form  of  a 

g 

hypergeometric  series  .  When  waveguide  fabrication  is 
done  by  diffusion,  surface  irregularities  can  be  thought 
of  as  causing  an  uneven  diffusion  depth.  This  will  yield 
an  index  profile  which,  under  sinusoidal  imperfections, 
takes  the  following  form: 

n(x,z)  =  ng  +An  cosh  [(x  -  xq  -  a  sin@z)/d]  (4) 

where  0  is  the  wave-vector  of  the  perturbation  and  is  given 
by  0  =  2ir/A  where  A  is  the  assumed  periodicity  of  the  perturba¬ 
tion.  It  is  assumed  that  the  ripple  amplitude  a  is  much 
smaller  than  the  effective  thickness  of  the  guiding  layer. 
Actual  radiation  modes  of  this  profile  are  complex  oscillatory 
hypergeometric  series.  However,  when  An<<ng,  the  exact 

solutions  approach  pure  oscillatory  solutions.  Approximating 
with  sinusoidal  fields  greatly  simplifies  calculations  without 
introducing  substantial  error.  This  assumption  has  been 
shown  valid  for  step-index  fibers  and  for  waveguides  with  small 
index  discontinuities9'10.  We  have  solved  u"der  single 
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mode  operation  just  below  the  cutoff  of  the  next  order  mode 
and  obtained  the  following  normalized  power  loss  in  the  range 
0  <  1  P  I  <  ko  ns  : 

AP/A3  \=  2t  n„  (An)2  pd  (1  +  (pd)2)2  (5) 

■plrv  ~r~  5  — ~2 .  t  ;~T' 

a  L  cosh  (pdir/2) 

where  A  is  the  free  space  wavelenath,  L  is  the  length  of  the 

2  2  2 

guide.  The  quantity  p  =  (k  n  )  -  (B  -  0)  where  8  is 

y  o  s  o  o 

the  propagation  constant  of  the  guided  mode  and  kQ  is  the 
propagation  constant  in  free-space. 

Next,  we  consider  diffused  surface  waveguides  with  an 
index  profile  given  by: 


n(x) 


n  + 
s 


in  e_<*/a> 


n 


air 


x  >,  t) 
x<  0 


(6) 


In  this  case,  imperfections  are  represented  by  sinusoidal 
perturbations  of  the  air-substrate  surface  such  that  n(x,z) 
of  the  imperfect  guide  is: 


n (x, z) 


n  +  An  e 
s 


- (x  +  a  sin  9 z) /d) 


(7) 


The  evaluation  of  the  losses  as  defined  in  (1)  was  carried  out 
by  using  the  exact  radiation  fields  rather  than  the  oscillatory 
approximation.  However,  it  was  assumed  that  the  radiation 
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fields  become  sinusoidal  as  x+°°.  Furthermore,  because  of 

the  large  air-dielectric  discontinuity,  we  consider  scattering 

only  into  the  substrate  as  the  predominant  type  of  radiation 

loss.  The  normal  radiation  loss  in  the  range 

k  n_  <  I  31  <  k  n  becomes: 
o  2  1  I  os 


AP  /  X 3  .  2“’  3<"s2  -  n22’2  Jv2<V)  lJi„  (V»  '  2 

P  a2L  uns 


lim 

T-+-0 


(|J-  <T)'2  +  fu+1(T1 


2-  T  Re  (iJ.  (T)  J*  ^ 

—  ip  ip  +  1 


where  V  = 


v  = 


b  = 


2d  k  J 2n  An 
o  v  s 

V  fh 


2  2 
n  Ci.-n 
ef  f  s 

2n  An 
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-  2d('Vs)2  -  2)1/2 


Equations  (5)  and  (8)  were  plotted  in  Fig.  1  illustrating 

the  radiation  losses  in  buried  and  diffused  surface  waveguides. 

In  order  to  numerically  evaluate  comparable  waveguides,  the 

same  propagation  vector  was  chosen  for  both  cases,  consistent 

with  single  mode  operation  just  below  the  cutoff  of  the  next 

order  mode.  The  parameters  were  n  =  1.5,  An  =  0.0035  and 

s 

kod  =  13 * 8  for  the  buried  waveguide,  and  ng  =  1.5,  An  =  0.01 


(T) 
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and  k  d  =  15.9  for  the  diffused  surface  waveauide. 
o 

It  can  be  seen  that  over  the  entire  range  of  periodici¬ 
ties  for  the  surface  roughness,  the  buried  guide  is  at 
least  one  order  of  magnitude  less  lossy  than  the  diffused 
guide.  Furthermore,  the  losses  of  the  buried  guide  are 
negligibly  small  for  short  range  roughness.  This  indicates 
that  burying  the  guiding  layer  results  in  a  significant 
reduction  in  scattering  losses  and  is  therefore  advantageous 

when  scattering  loss  is  of  major  concern.  This  is  in  agree- 

2  3 

ment  with  the  experimental  results  in  glass  ’ 

For  completeness,  we  compare  these  losses  to  those  of 

abrupt  uniform  waveguides,  symmetric  and  asymmetric,  again 

requiring  the  same  propagation  vector  in  each  waveguide. 

5-7 

These  losses  were  calculated  from  equations  of  Marcuse 
In  order  to  make  correct  analogy  with  diffused  and  buried 
guides,  we  assumed  both  boundaries  of  the  abrupt  guide  were 
perturbed  in  parallel  with  the  same  periodicity.  This  cor¬ 
responds  to  a  =  0  in  Marcuse's  notation,"*  and  has  been  termed 
microbending"*.  These  results  are  shown  as  the  solid  lines  in 
Fig.  2.  In  order  to  make  numerical  comparison  with  the 
results  in  Fig.  1,  the  same  substrate  refractive  index  propaga¬ 
tion  constant  were  assumed.  By  comparison  with  Fig.  1,  it  can 
be  seen  that  the  peak  loss  of  the  symmetric  abrupt  guide  is 
similar  to  that  of  the  buried  guide,  and  the  peak  loss  of  the 
asymmetric  abrupt  guide  is  also  similar  to  that  of  the  surface 
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lized  radiation  power  loss  of  equivalent  surface  diffused  and  buried 
uides  wit-b  sinusoidal  surface  irregularity.  Included  on  the  figure 
ketches  of  the  refractive  index  profiles. 


A  A 


^/V  O'OOL 


diffused  guide.  In  general,  it  can  be  said  that  abrupt  guides 
have  comparable  losses  to  their  equivalent  graded  index  guides. 

The  abrupt  guides,  however,  can  be  seen  to  have  higher 
losses  at  particular  values  of  A  corresponding  to  resonances 
of  the  leaky  modes  of  the  radiation  fields.  $uch  leaky  modes 
do  not  exist  in  graded  index  waveguides,  which  are  therefore 
less  sensitive  to  short  wavelength  perturbations.  This  result 
was  also  observed  by  Marcuse  in  microbend  loss  calculations 
of  slab  guides  with  a  parabolic  index  profile^.  This  reduc¬ 
tion  of  radiation  loss  at  small  periodicities  of  surface 
roughness  in  graded  index  guides  is  significant  since  short 

range  random  roughness  is  the  predominant  type  on  most 

12 

polished  surfaces 

It  should  be  mentioned  that  for  abrupt  waveguides  fab¬ 
ricated  by  various  deposition  techniques,  there  is  a  good 
possibility  that  the  perturbation  of  the  two  walls  is  uncor¬ 
related.  To  show  the  effect  of  a  phase  difference  between 
the  perturbation  in  both  walls,  the  radiation  loss  for  the 
symmetric  case  in  which  the  waveguide  walls  are  perturbed 
symmetrically  about  the  center  of  the  waveguide  is  plotted  as 
the  dotted  line  in  Fig.  2.  It  can  be  seen  that  in  this  case 
the  guide  is  more  sensitive  to  long  period  wall  roughness 
than  when  the  perturbations  occur  in  parallel.  Uncorrelated 
surfaces  will,  in  general,  have  losses  somewhere  between  the 
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two  extremes.  The  situation  when  the  two  boundaries  are 
uncorrelated  and  the  loss  is  correspondingly  higher  is  more 
likely  to  be  found  in  abrupt  films  than  in  graded  guides, 
making  graded  layers  more  attractive  in  this  respect. 

In  conclusion,  we  have  shown  quantitatively  that  buried 
graded  index  guides  have,  in  principle,  lower  losses  due  to 
surface  roughness  than  do  surface  diffused.  It  is  therefore 
clear  that  techniques  must  be  developed  to  bury  the  guiding 
layer  for  low  loss  integrated  optic  circuitry. 
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c .  Absorption  Loss  in  LiNbO^  Waveguides 

Absorption  loss  was  measured  in  LiNbO^  waveguides  by  using 
calorimetry.  This  technique  was  demonstrated  to  be  excellent 
and  very  sensitive  method  for  determining  this  loss  mechanism. 
The  results  of  our  first  work  was  published  in  Applied  Physics 
Letters.  A  reprint  of  this  work  is  included  in  this  section, 
describing  the  technique  and  tabulating  the  losses  measured  in 
LiNbO^  outdif fused  waveguides.  Later  measurements  were  made  in 
LiNb03  Ti  in-diffused  waveguides  as  well  as  with  another  series 
in  Table  I.  More  complete  details  on  comparison  of  mechanisms 
for  loss  is  contained  in  Section  d. 

TABLE  3-1 

Absorption  Loss  in  LiNbO^  Waveguides 
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The  results  of  these  measurements  can  be  summarized  by 

the  following  facts: 

1.  The  absorption  loss  in  waveguides  is  much  smaller  than  the 
scattering  loss,  for  both  in-diffused  and  out-diftused 
waveguides . 

2.  The  TM  mode  has  more  absorption  than  the  TE  mode  (in- 
diffused  waveguides).  The  reason  for  this  is  not  understood 
at  this  time.  Nevertheless,  these  losses  are  much  smaller 
than  the  total  loss,  and  at  the  present  time  are  negligible. 

3.  The  TE  mode  has  the  same  loss  whether  the  waveguide  is  fabric¬ 
ated  by  out-diffusion  or  by  in-diffusion  of  Titanium.  This 
means  that  the  Ti  ions  do  not  introduce  significant  absorption. 

4.  The  absorption  loss  in  TE  modes  is  very  close  to  that  in 
bulk  LiNbO-j.  Because  of  the  errors  in  the  measurements 
and  the  small  value  of  the  absorption  loss,  it  is  not  known 
whether  the  loss  in  waveguides  is  the  same  as  in  bulk,  or 
bigger  by  a  factor  of  two.  Nonetheless,  the  absorption 
losses  are  an  order  of  magnitude  smaller  than  the  total  loss. 

5.  The  loss  at  1.06  ym  is  three  times  smaller  than  at  0.63  ym. 

6.  There  is  a  dichroism  in  bulk  LiNb03,  which  exhibits  twice  the 
absorption  for  the  extraordinary  polarization  than  for  the 
ordinary  polarization.  This  effect  has  'not  been  reported 
before.  Experiments  are  underway  to  study  the  dichroism  at 
other  wavelengths. 

7.  The  total  loss  in  LiNbO^  bulk  is  about  twice  the  absorption 
loss.  The  difference  being  presumably  due  to  scattering. 


8.  The  total  loss  in  LiNbO^  waveguides  is  about  an  order  of 
magnitude  bigger  than  the  total  loss  in  bulk  LiNbO^- 
This  is  the  first  time  that  several  of  these  facts  have 
been  reported,  and  not  all  these  properties  can  be  explained 
at  the  present  time.  However,  clearly  absorption  loss  is  not 
a  dominant  mechanism  in  LiNbO^  waveguides. 
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Calorimetric  measurement  of  LiNb03  waveguide  absorption 
losses 
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This  paper  reports  the  first  measurement  of  optical-absorption  loss  in  LiNbO, 
waveguides  and  bulk  material.  This  is  also  the  first  demonstration  of  calorimetry  as  a 
feasible  technique  for  measuring  waveguide  absorption  in  integrated  optics.  Results 
indicate  that  in  LiNbO,  out-diffused  waveguides  the  absorption  is  two  orders  of 
magnitude  lower  than  the  total  loss,  implying  that  the  predominant  loss  mechanism  is 
scattering.  The  ultimate  sensitivity  limit  for  calorimetric  determination  of  absorption 
losses  appears  to  be  —  10  *  dB /cm  for  optimized  geometries. 

PACS  numbers.  42.80.Lt,  42.80.Sa,  42.82.  +  n,  42.70.Fh 


Measurement  of  absorption  losses  by  laser  calorimetric 
techmques  is  common  in  the  fields  of  high-power  laser  op¬ 
tics1  and  fiber  optics.1  The  application  of  this  technique  to 
determine  the  absorption  loss  in  LiNbO,  waveguides  and 
bulk  material  at  1.06 pm  is  reported.  The  results  are  com¬ 
pared  to  total  loss  measurements  at  this  wavelength  in  order 
to  elucidate  the  loss  mechanisms. 

In  a  calorimetric  measurement  the  temperature  rise  of  a 
sample  during  laser  irradiation  is  detected.  Such  data  per¬ 
mits  a  sensitive  and  absolute  determination  of  the  absorption 
loss  in  a  sample.  Bulk  absorption  coefficients  of  less  than 
2  -  10  ’  cm  ( 10  '  dB/cm)  have  been  measured  at  infrared 
wavelengths.1  In  contrast  to  transmission  measurements,  ca- 
lorimelrically  determined  waveguide  losses  depend  only  to 
second  order  on  insertion  losses.  Since  the  absorption  is  inde¬ 
pendent  of  the  laser-beam  diameter  and  since  a  waveguide 
will  confine  the  beam  along  the  entire  length  of  the  sample, 
the  sensitivity  of  a  waveguide-loss  measurement  can  be  com¬ 
parable  to  that  of  a  bulk-loss  measurement  in  samples  of 
equal  thermal  mass.  Measurements  in  samples  of  different 
geometries  are  reported,  and  we  estimate  that  for  proper 
choice  of  sample  geometry,  the  minimum  measurable  wave¬ 
guide  absorptivity  is  only  10  times  larger  than  the  minimum 
measurable  bulk  absorptivity  We  therefore  conclude  that 
laser  calorimetry  is  capable  of  measuring  much  lower  ab¬ 
sorption  losses  than  any  of  the  integrated-optics  techniques 
in  current  use. 

LiNbO,  is  an  important  material  for  active  devices  in 
integrated  optics  and  is  known  to  form  low-loss  waveguides 
by  out-diffusion  and  Ti  in-diffusion.* '  Much  of  the  work  to 
date  on  losses  in  LiNbO,  waveguides  has  been  directed  to¬ 
ward  understanding  optical  damage  effects  in  the  visible.'  At 
1.06/tm  there  is  negligible  optical  damage,  and  this  wave¬ 
length  region  is  of  particular  interest  because  of  the  mini¬ 
mum  optical-fiber  loss.  Total  losses  (absorption  plus  scatter¬ 
ing)  in  LiNbO,  bulk  material  have  been  shown  to  be  less  than 
0. 1  dB/cm  at  1. 1 5/rm.’ Total  losses  in  a  Ti  in-diffused  wave¬ 
guide  have  been  measured  at  0.63  /cm  (at  power  levels  below 
the  optical  damage  threshold)  and  yielded  a  value  of  1 .0 
dB/cm.’  In  this  paper,  we  report  the  losses  at  l  .Obyrm  in 


waveguides  prepared  on  y-cut  LiNbO,  substrates  by  out-dif¬ 
fusion  in  flowing  oxygen  for  1  h  at  950  “C. 

The  laser  used  in  this  experiment  was  Nd  :  YAG  oper¬ 
ated  cw  at  1.06/jm  and  producing  200  mW  of  TEM«  inci¬ 
dent  power.  No  optical  damage  effects  were  observed  at  this 
power  level  The  initial  laser-beam  diameter  was  0.76  mm 
FWHM  and  the  polarization  was  aligned  to  excite  the  TE 
mode  of  the  waveguide.  The  focusing  optics  consisted  of  a 
cylindrical  lens  (f.l.  -  80  mm)  and  a  microscope  objective 
(10X  ,  N.A.  =  0.25).  A  simple  lens  (f.l.  =  40  mm)  and  a  cy¬ 
lindrical  lens  (f.l  =  40  mm)  collected  the  output  of  the 
waveguide  and  focused  it  onto  a  power  meter  with  a  slit 
positioned  to  block  scattered  light.  For  measurements  of 
bulk  absorption  and  transmission,  the  slit  was  removed.  The 
samples  were  mounted  on  glass  capillary  tubing  with  cyan- 
oacrylic  cement  to  minimize  conduction  losses.  Small-diam¬ 
eter  copper/Constantan  thermocouples  (0.003  in..  Omega 
Engineering)  were  similarly  glued  to  the  samples.  The  refer¬ 
ence  thermocouple  was  embedded  in  an  aluminum  block.  A 
Plexiglas  enclosure  provided  isolation  of  the  sample  and 
thermocouples  from  air  currents.  Relative  sample  tempera¬ 
ture  changes  before,  during,  and  after  laser  irradiation  were 
monitored  with  a  nanovoltmeter  (Keithley  147). 

In  order  to  explore  the  usefulness  of  calorimetry  as  a 
measurement  technique  for  integrated  optics,  we  measured 
the  losses  in  two  very  different  geometries,  representing  the 
worst  and  best  cases  The  results  are  summarized  in  Table  I. 
Sample  A  (dimensions  0.3  x  2.4  X  2.5  cm,  m  -  8.6  g)  had  a 
relatively  large  cross-sectional  area  and  a  long  thermal  time 
constant,  which  produced  a  heating  curve  linear  with  time, 
as  shown  in  Fig.  1(a).  The  absorption  coefficient  was  deter¬ 
mined  as  follows*: 


where  m  is  the  mass  of  the  sample,  C  is  the  specific  heat,  /  is 
the  length  of  the  sample,  and  P„  is  the  incident  laser  power. 
The  quantities  T  (  and  T,  are,  respectively,  the  slopes  of  the 
temperature-vs-time  curves  and  the  temperatures  of  the 
samples  for  the  three  regions  shown  in  Fig.  1(a). 
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FIG  1  Temperature-vs-time  bahavtor  for  the  two  sample  geometries. 


Sample  B  (dimensions  0. 1  x  0. 1  x  2. 5  cm,  m  =  0. 1 4  g), 
on  the  other  hand,  had  a  small  cross-sectional  area  and  a 
thermal  behavior  which  was  dominated  by  convection  and 
radiation  losses,  as  shown  in  Fig  1(b).  For  this  case,  the 
absorption  coefficient  was  given  by 

a  (cm  ')  =  mCAT /lrPt),  (2) 

jvhered  Tis  the  steady-state  sample  temperature  rise  and  r  is 
the  measured  thermal  time  constant  of  the  sample.”  The  val¬ 
ues  of  the  absorption  coefficients  obtained  for  the  two  differ¬ 
ent  geometries  agree  within  experimental  error,  and  the 
waveguide  loss  is  only  about  twice  that  observed  in  the  bulk 
This  indicates  that  out-diffusion  does  not  introduce  signifi¬ 
cant  absorption  centers  in  LiNbO,  at  1  06/zm.  Both  the  bulk 
and  waveguide  absorption  losses  are  considerably  lower 
than  the  total  losses,  as  will  be  discussed  below. 


The  total  loss  in  the  optical  system  through  either  the 
bulk  or  the  waveguide  can  be  described  by 


jj_  _  Lf  (  (1  -  R  expt  -  («  +  P  V  1  V 
Pn  '  \  1  —  R  2  exp[  —  2(a  +  0)1  ]  ' 


(3) 


where  F,  and  P0  are  the  transmitted  and  incident  powers, 
respectively,  L0  and  L ,  are  the  fractional  transmissions  of 
the  focusing  and  collecting  optics,  respectively,  is  the 
waveguide  coupling  efficiency,  R  is  the  reflectivity  for  a  sin¬ 
gle  interface,  a  and  p  are  the  length-dependent  absorption 
and  scattering  losses,  respectively,  and  /  is  the  sample  length. 
The  term  in  large  parentheses  in  Eq.  (3)  is  the  equation  for 
loss  in  a  sample  including  noncoherent  multiple  reflec¬ 
tions.  10  The  transmission  coefficients  of  the  input  and  output 
lenses  were  L0  =  0.73  and  Lr  =  0.81 ,  The  coupling  efficien¬ 
cy  was  determined  by  calculating  the  overlap  between  the 
10-/im-deep  waveguide  and  the  incident  spot  size  (9.8  X  23 
^im).  The  coupling  efficiency  was  >0.90  with  the  lower 
bound  used  in  the  data  reduction. 


The  total  loss  in  bulk  LiNbO,  (0.13  dB/cm,  sample  A) 
is  comparable  to  that  reported  by  Kaminow  and  Stulz’  at 
1.1 5  urn  (0.1  dB/cm).  In  both  cases,  this  is  only  an  upper 
bound  on  the  true  total  bulk  loss.  The  total  waveguide  loss  is 
approximately  the  same  for  both  samples  and  almost  two 
orders  of  magnitude  greater  than  the  absorption  loss.  This 
indicates  that  the  dominant  optical-loss  mechanism  in  these 
waveguides  is  scattering.  It  is  interesting  to  note  that  the 
total  waveguide  loss  at  1.06  gm  is  comparable  to  the  loss 


reported  for  a  Ti  in-diffused  waveguide  at  0 1>3  gm  using  an 
interferometric  technique  (1.0  dB/cm). 

For  the  optimized  geometry  of  small  rodlike  samples 
whose  length  is  much  greater  than  both  the  width  and  the 
height,  the  minimum  measurable  absorptivity  is  proportion¬ 
al  to  the  perimeter  of  the  entrance  face  Therefore,  only  small 
gains  in  sensitivity  are  made  by  decreasing  the  size  of  the 
sample.  The  minimum  measurable  absorption  coefficient  of 
samples  similar  in  size  to  sample  B  is  estimated  to  be  5  \  10  4 
dB/cm.  The  sensitivity  in  this  ease  w  ill  be  limited  by  scatter¬ 
ing  directly  onto  the  thermocouple  and  may  be  further  im¬ 
proved  by  optimizing  the  coupling  efficiency  of  the  laser  into 
the  waveguide.  For  large  samples  whose  thermal  behavior  is 
best  analyzed  by  Eq.  (1),  the  minimum  measurable  absorp¬ 
tivity  is  proportional  to  the  cross-sectional  area  of  the  en¬ 
trance  face  of  the  substrate.  As  a  rule  of  thumb,  for  experi¬ 
mental  condition  similar  to  those  in  this  work,  the  minimum 
measurable  absorptivity  is 

amm  (cm  ')  =  3x  10  A. 

where  A  is  the  area  in  cm/  In  generai,  the  sensitivity  of  the 
calorimetric  measurements  could  be  improved  for  both  the 
large  and  small  samples  by  (a)  increasing  the  thermal  isola¬ 
tion  by  placing  the  samples  in  a  vacuum;  (b)  improving  the 
laser  stability;  (c)  eliminating  unguided  light. 

In  conclusion,  we  have  demonstrated  the  usefulness  of 
the  calorimetric  technique  for  measuring  absorption  loss  in 
planar  waveguides.  It  is  estimated  that  the  minimum  wave¬ 
guide  absorption  coefficient  measurable  for  samples  with  op¬ 
timum  geometry  is  on  the  order  of  10  4  dB/cm.  For  the  spe¬ 
cific  example  of  LiNbO,  out-diffused  waveguides,  the 
absorption  coefficient  at  1.06yum  is  2  x  10  1  dB/cm.  This 
result  is  to  be  compared  with  the  first  measurement  of  ab¬ 
sorption  coefficient  in  bulk  LiNbO,,  1  x  10  1  dB/cm  report¬ 
ed  here.  The  total  losses  measured  in  both  the  bulk  and  wave¬ 
guide  were  found  to  be  comparable  with  previously  reported 
results.’  These  results  show  that  in  the  waveguide  the  ab¬ 
sorption  loss  is  considerably  less  than  the  total  loss,  giving 
the  first  indication  that  the  absorption  losses  in  LiNbO,  are 
negligible  at  !.06^m  and  demonstrating  that  scattering  is 
the  limiting  mechanism  in  the  waveguide  transmission. 
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d.  Loss  Mechanisms  in  LiNbO^  Waveguides 

Loss  mechanisms  in  optical  waveguides  can  be  classified 
into  three  types;  absorption,  volume  scattering  due  to  refrac¬ 
tive  index  irregularities  within  the  interior  of  the  guide, 
and  surface  scattering  due  to  waveguide  boundary  imperfections. 

We  have  separated  these  mechanisms  and  proved  that  the  latter 
is  the  main  source  of  length-dependent  loss  in  LiNbO^  waveguides. 

We  showed  in  Section  c  that  absorption  is  negligibly  small 
compared  to  the  overall  guide  loss.  It  remains,  therefore,  to 
differentiate  between  volume  and  surface  scattering  loss  mech¬ 
anisms  for  length  dependent  losses. 

The  basic  technique  used  for  differentiating  between  scat¬ 
tering  due  to  surface  roughness  and  scattering  due  to  volume 
inhomogenienities  is  the  comparison  of  in-plane  and  out-of-plane 
scatter.  Scattering  from  volume  inhomogenieties  will  be  roughly 
isotropic,  while  scattering  from  surface  roughness  will  cause 
primarily  out-of-plane  scatter.  This  latter  face  is  because, 
for  this  case,  one-dimensional  phase  matching  occurs,  in  which 
surface  roughness  couples  guided  modes  into  small-angle  radiation 
modes. 
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These  modes  radiate  out-of-plane  and  are  much  stronger  than  the 
in-plane  radiation  caused  by  surface  roughness. 

Estimations  of  in-plane  scatter  are  relatively  simple  to 
obtain,  because  the  beam  profile  of  light  travelling  in  the 
waveguide  can  be  compared  to  the  beam  profile  of  light 
travelling  in  the  bulk  of  the  sample.  Any  broadening,  or  wings 
on  the  profile  of  the  guided  light  will  be  due  to  scatter 
within  the  waveguide.  In  order  to  obtain  these  profiles,  LiNbO^ 
samples  were  polished  and  end-fire  coupled.  The  near  field  of 
the  exit  face  was  imaged  onto  a  scanning  detector.  Intensity  scans 
were  taken  of  the  light  when  it  was  guided  and  compared  to  those 
when  the  light  was  traveling  through  the  bulk  substrate. 

The  overlap  integral  between  the  two  profiles  determined 
the  fraction  of  in-plane  scattering  of  the  guide  relative  to 
that  of  the  bulk.  Measurements  of  LiNbC>3  out-diffused  guides 
showed  that  the  in-plane  scatter  was  less  than  0.05  dB/cm, 
proving  that  volume  scatter  is  not  a  predominant  loss  mechanism 
in  these  waveguides. 

All  LiNbO^  waveguides  were  studied  in  with  end-fire 
coupling,  in  the  configuration  with  which  they  will  be  used  for 
fiber  optics  switches.  Typical  near-field  profiles  of  the  exit 
light  are  shown  in  Fig.  3-1.  It  can  be  seen  that  there  is  a 
considerable  amount  of  light  in  both  the  substrate  and  the  air 
above  the  waveguide.  The  ratio  of  the  energy  inside  the  wave¬ 
guide  to  that  in  the  wings  gives  a  measure  of  the  insertion  loss 
of  the  waveguide.  This  insertion  loss  will  be  due  to  the 


mode-mismatch  between  the  incident  light  and  the  waveguide  and 
also  due  to  small  angle  scatter  loss  in  the  waveguide.  Typical 
numbers  for  the  fraction  of  light  in  each  regime  is  shown  in 
the  first  column  of  Table  II,  for  both  out-diffused  and 
in-diffused  waveguides.  Included  in  the  table  is  the  amount  of 
light  which  was  undetected.  All  measurements  were  normalized 
to  the  transmission  in  the  bulk.  Since  bulk  losses  were  an 
order  of  magnitude  smaller  than  waveguide  losses  (section  c) , 
and  we  wanted  to  eliminate  the  need  to  calculate  reflection 
losses,  the  normalization  to  the  bulk  transmission  gave  a 
simple  way  to  do  this.  It  can  be  seen  that  the  undetected 
light  (lost  due  to  absorption  or  large-angle  scatter)  is  very 
much  smaller  than  small  angle  scatter  or  input  coupling  losses. 

It  was  possible  to  get  a  direct  measure  of  the  light  which 
was  scattered  into  the  air  by  surface  roughness.  We  placed  a 
block  above  the  guide  at  the  entrance  face  to  eliminate  all 
light  due  to  mode  mismatch,  and  measured  the  output  in  the  air 
due  to  surface  scatter.  The  difference  between  this  signal  and 
the  signal  measured  without  the  block  represents  a  measure  of 
the  coupling  mismatch  or  small  angle  scattering  into  the  air 
due  to  roughness  of  the  entrance  edge.  In  order  to  estimate 
the  difference  between  these  two  quantities,  we  assumed  a 
Gaussian  mode  input  and  a  Bessel  function  guided  wave  and 
calculated  the  overlap.  The  results  are  listed  in  the  table 


TABLE  II 


Contributions  of  Various  Mechanisms  at  the  Exit  Face  of 
LiNbO^  Waveguides  1  inch  Long  at  .633  pm.  Measurements  normalized  to 

bulk  transmission 


OUT-DIFFUSED 

Total 

Measured 

/t  =  10 
yTE  mode 

um  Guided 

5  pm 

mode  width 

In  Guide 

43.5% 

Small  Angle 

Scatter 

Coupling 

Mismatch 

Entrance 

Edge  loss 

In  Air 

9% 

3% 

(3%) 

(3%) 

In  Substrate 

46% 

43% 

(4%) 

(<3%) 

Undetected 

Light 

1.5% 

IN-DIFFUSED 

/  t  =  3  pm 
y  TE  mode 

Guided 

mode  width 

2.6  ym 

In  Guide 

33% 

In  Air 

14% 

6.5% 

(3.5%) 

(4%) 

In  Substrate 

50% 

43% 

(7%) 

(<  3% ) 

Undetected 

Light 

3% 

under  "coupling  mismatch"  and  are  placed  in  brackets  to  show 
that  they  are  the  result  of  calculations .  The  light  which  is 
in  the  air  and  not  represented  by  the  coupled  mode  mismatch 
or  by  the  small  angle  scatter,  is  due  to  an  entrance  edge 
scatter,  and  is  listed  in  the  last  column. 

It  is  not  possible  to  make  separate  measurements  in  the 
substrate.  However,  we  made  a  few  gross  assumptions  to  develop 
some  estimates.  We  calculated  the  light  in  the  substrate  due 
to  the  expected  mode  mismatch,  and  then  assumed  that  the  rest 
of  the  light  in  the  substrate  was  due  to  small  angle  scatter. 
These  numbers  are  included  in  the  table,  in  brackets.  Also 
included  is  the  fact  that  we  assumed  the  entrance  edge 
scatter  into  the  substrate  to  be  less  than  that  in  the  air. 

With  the  assumption  that  the  coupling  mismatch  and  entrance 
edge  loss  all  occur  at  the  beginning  of  the  waveguide,  before 
length/dependent  scatter  loss  becomes  appreciable,  the  length- 
dependent  small  angle  scatter  loss  can  be  calculated.  For 
out-diffused  waveguides,  the  initial  transmission  at  the  front 
edge  was  87%,  corresponding  to  a  0.6  dB  loss  (excluding 
reflections).  The  remainder  of  the  total  3.6  dB  loss  is  due 
to  small  angle  scattering,  yielding  a  length-dependent  loss 
of  1.1  dB/cm.  For  the  in-diffused  guide,  the  transmission  loss 
of  the  front  edge  is  17.5%  corresponding  to  a  0.83  dB  loss. 

The  remainder  of  the  total  4.8  dB  loss  is  distributed  along  the 
length  and  corresponds  to  a  loss  of  1.6  dB/cm.  Finally,  the 


reflectivity  has  not  been  included  in  this  analysis,  since 
these  measurements  were  all  made  relative  to  the  bulk  trans¬ 
mission.  The  reflection  losses  for  two  faces  of  LiNb03, 
with  a  reflectivity  of  15%  each  is  1.3  dB.  This  means  that 
the  total  insertion  loss  for  the  two  samples  was  5  and  6  dB 
respectively.  The  samples  were  carefully  prepared  to  have  the 
same  waveguide  thickness  and  good  quality  edges,  so  that 
rounding  and  chips  were  eliminated.  The  plot  of  these  losses 
as  a  function  of  waveguide  length  yields  the  loss  per  unit 
length.  The  results  are  shown  in  Fig.  3-2,  for  samples  up  to 
7.5  cm  in  lenqth. 

The  results  of  the  longest  sample  do  not  belong  on  the  same 
graph  since  that  guided  mode  was  a  different  thickness.  We 
have  included  the  measurements  for  completeness,  however. 

The  best  estimate  we  have  for  the  length  dependence  of  the  loss 

is  drawn  in  as  a  straight  line  and  corresponds  to  a  =  0.6  dB/cm. 
The  results  are  also  consistent  with  a  =  0.3  dB/cm.  These 

numbers  are  somewhat  lower  than  obtained  previously.  The 
difference  is  real,  because  these  waveguides  were  thicker  than 
those  studied  previously  and  will  therefore  have  less  loss. 

The  insertion  loss  at  zero  thickness  includes  the  Fresnel 
reflection  loss  of  1.3  dB.  The  remainder  of  3.5  dB  is  due  to 
any  mode-mismatch  and  edge-coupling  loss. 

The  results  of  this  study  have  indicated  that  the  mech¬ 
anisms  for  loss  in  end-coupled  LiNbO^  waveguides  are  both 
length-dependent  losses  due  to  scatter  from  the  waveguide  surface 
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qure  3-2.  Insertion  loss  of  out-diffused  LiNbC^  waveguides  of  different  length 
Samples  marked  •  were  pure  out-diffused.  Samples  marked  X  weie  meas 
the  presence  of  in-diffused  titanium.  Guided  modes  were  10  pm  thick 
15 1 —  Samples  marked  a  had  guided  mode  profiles  5  urn  thick. 


and  also  losses  due  to  poor  coupling,  mode  mismatch  and 
reflection  at  the  entrance  edge.  The  first  set  of  measure¬ 
ments  indicated  a  relatively  large  value  for  the  length- 
dependent  loss,  making  this  the  predominant  loss  mechanism. 

The  second  set  of  measurements  indicated  that  the  edge  losses 
were,  in  general,  larger  than  the  length  dependent  losses. 

More  work  is  needed  to  sort  out  these  effects  in  order  to 
reduce  the  loss. 

e.  Techniques  to  Reduce  LiNbO-,  Waveguide  Insertion  Loss 

In  order  to  determine  whether  the  primary  source  of  loss 
was  due  to  surface  scatter  and  to  reduce  this  effect  we  tried 
several  experiments.  First  we  tried  polishing  the  top  of  the 
samples.  It  is  well  known  that  in-diffusion  of  Titanium 
causes  the  surface  of  the  sample  to  be  raised  slightly,  and 
a  surface  roughness  to  appear.  We  therefore  decided  to 
measure  the  losses  before  and  after  the  surface  was  polished. 
The  results  are  shown  in  Table  III.  It  can  be  seen  that  the 
waveguide  became  10  ym  thinner  during  the  polishing  and  the 
waveguide  losses  typically  increased.  Unambiguous  inter¬ 
pretation  of  these  measurements  is  not  possible  since  there 
were  several  mechanisms  which  may  have  contributed  to  higher 
losses  in  the  thinner  waveguides  after  polish.  First,  thinner 
waveguides  typically  have  higher  scatter  losses.  Second,  thin¬ 
ner  waveguides  have  higher  input  coupling  losses.  Both  these 
effects  were  reported  in  Section  d. 


The  losses  of  the  samples  reported  in  table  ITT  ire 
hiqher  than  those  reported  in  section  d.  These  wo- e  some  of 
the  first  measurements  we  made,  and  it  was  discovered  after 
these  measurements  that  the  corners  of  the  waveguides  tended 
to  be  rounded.  Commercially  polished  edges  were  found  to  be 
inadequate  for  good  input  coupling  efficiencies.  Careful 
polishing  procedures  were  introduced  on  the  front  edge  of 
the  samples  and  the  insertion  losses  were  reduced  to  typically 
3  dB . 

In  order  to  establish  definitively  the  success  of  sur¬ 
face  polish,  new  measurements  should  be  made  with  samples 
whose  edges  have  been  properly  prepared  and  also  which  have 
not  been  thinned  to  such  an  extent.  It  was  not  possible  to 
make  these  measurements  before  the  end  of  the  contract. 

In  order  to  further  investigate  the  cause  of  the  loss  in 
samples  with  good  end  faces  another  experiment  was  performed. 
Losses  due  to  surface  scatter  depend  on  the  index  difference 
between  the  waveguide  and  the  medium  above  the  surface.  If 
the  medium  above  the  surface  is  changed  from  air  to  a  trans¬ 
parent  medium  with  a  higher  index,  the  amount  of  surface 
scatter  should  be  reduced.  We  performed  this  experiment  by 
measuring  the  losses  before  and  after  coating  the  top  of  a 
I.iNbO-j  waveguide  with  a  polymer. 

The  results  are  shown  in  Table  3-IV.  The  polymer  was  in 
‘  * •  form  of  a  glue  with  refractive  index  1.4,  and  it  was 
i:  '  ■  ir-d  carefully  in  order  not  to  disturb  the  entrance  and 
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TABLE  3-III 


Insertion  loss  (reflection  loss  included)  of  TE  modes  at 
0.63  pm  in  2.54  cm  LiNb03  waveguide  before  and  after  10 
minutes  of  hand  polishing  with  Syton  on  a  silk  pad. 


Before  After 


Sample 

Number 

Guided  Mode 
Thickness 
Micrometers 

Loss 

dB 

Guided  Mode 
Thickness 
Micrometers 

Loss 

dB 

A 

20 

9.0 

10 

15.0 

B 

25 

9.0 

10 

8.8 

C 

30 

7.9 

30 

9.4 

D 

30 

5.0 

20 

6.2 

Insertion  Loss  Change  of  TE  Guided  Mode  for  Y-cut  Propagating  LiNb03  Waveguides 
The  cladding  refractive  index  was  1  without  glue  and  the  glue  index  was  1.4 
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exit  faces  to  the  waveguide,  lest  it  affect  the  input  and  output 
coupling  efficiencies.  All  the  measurements  listed  in  Table  3-TV 
were  consistent  with  no  change  in  transmission  due  to  the  polymer 
in  the  waveguide  surface.  The  error  for  these  measurements  is 
about  10%,  expressed  as  a  percentage  of  the  length-dependent 
loss  (assumed  to  be  0.6  dB/cm) .  This  number  can  be  compared  to 
that  expected  from  theory. 

It  was  shown  in  Section  b  (Eq.  8)  that  the  normalized  radia¬ 


tion  loss  due  to  surface  roughness  is 


0  3  .  2  2. 

2r  (ns  -  nc  ) 


where  F  is  a  function  which  depends  on  the  details  of  the  mode  profile 

and  the  amplitude  and  periodicity  of  the  roughness.  The  quantity  ng 

is  the  index  of  the  substrate  and  nc  is  the  index  of  the  cladding. 

If  the  refractive  difference  between  the  guide  and  the  substrate 

is  much  smaller  than  the  difference  between  the  guide  and  air 

(or  polymer) ,  the  mode  profile  is  determined  primarily  by  the 

diffusion  profile  and  does  not  depend  on  the  index  of  the  cladding 

layer.  Under  this  approximation,  the  scatter  loss  should  be 

2 

2  2 

reduced  due  to  the  reduction  in  quantity  (n_  -  n^  )  .  This 

s  c 

difference  between  n  =  1  =  air  and  nc  =  polymer  is  40%.  We 
measured  -'10%  change. 

Our  inability  to  measure  a  40%  reduction  ir.  waveguide  loss 
may  be  due  to  several  causes.  First,  the  large  amount  of  inser¬ 
tion  loss  error  due  to  mode  mismatch,  reflections,  and  edge  defects 
may  have  swamped  the  length-dependent  losses  which  we  seek  to 
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reduce.  Second,  the  assumption  that  the  node  profile  does  not 
change  due  to  introduction  of  a  larger  index  cladding  layer  may 
be  inaccurate.  In  addition  to  making  the  theoretical  estimate 
wrong,  a  mode-change  due  to  the  polymer  cladding  will  cause  a 
mode-mismatch  between  light  initially  guided  in  an  air-cladded 
layer  and  the  mode  when  it  enters  the  polymer-clad  region.  The 
mode  mismatch  at  this  transition  may  introduce  additonal  loss 
masking  the  loss  reduction  we  intended  to  measure.  Thirdly, 
the  polymer  we  used  may  have  been  lossy  or  not  have  formed  a 
perfect  interface  between  the  waveguide  and  the  polymer.  Bubbles 
or  inhomogeneities  in  the  polymer  would  introduce  additional 
sources  of  loss.  Finally,  our  assumption  that  the  length- 
dependent  losses  are  due  to  scatter  from  the  surface  imperfec¬ 
tions  may  be  erroneous.  More  measurements  are  required  to  determine 
whether  or  not  reducing  the  index  discontinuity  between  the 
waveguide  and  cladding  layer  reduces  the  waveguide  loss. 

In  conclusion,  the  only  technique  we  have  found  which  sig¬ 
nificantly  reduces  the  insertion  loss  of  end-coupled  LiNbO^ 
waveguides  is  to  use  great  care  in  polishing  the  endfaces  of  the 
waveguides,  to  eliminate  any  rounding  of  the  corners.  Techniques 
attempted  to  reduce  the  loss  due  to  surface  scatter  did  not 
produce  any  noticeable  improvement. 

Clearly,  if  the  waveguide  could  be  buried  below  the  surface, 
the  requirements  on  the  polish  of  the  end  faces  would  be  reduced. 

In  addition,  any  loss  due  to  surface  roughness  would  be  reduced. 
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CHAPTER  4  .  STUDIES  OF  MULTIMODE  WAVEGUIDE  FABRICATION  IN 
C-PROPAGATING  LiNb03 

In  the  first  chapter  it  was  described  how  the  90°  distri¬ 
buted  Bragg  deflector  can  act  as  a  polarization  analyzer, 
deflecting  one  polarization  and  transmitting  the  other.  This 
is  the  integrated  optical  analog  of  the  Glan-Thompson  prism. 

This  suggests  its  use,  along  with  the  electro-optic  effect,  to 
make  a  multimode  integrated  optics  fiber  switch.  This  occurs 
when  phase  modulation  is  used  to  cause  polarization  modulation, 
converting  TE  to  TM  modes.  The  distributed  Bragg  deflector 
then  transmits  the  TE  polarization  and  deflects  the  TM  polari¬ 
zation.  If  two  exit  fibers  are  placed  as  shown  in  Fig.  4-1, 
an  integrated  optics  multimode  fiber  optics  switch  is  the  result 

Our  studies  of  the  beam  expander,  therefore,  led  to  an 
interest  in  demonstrating  the  device  in  LiNbC>3,in  conjunction 
with  an  integrated  optical  m  .ltimode  switch.  This,  in  turn, 
led  to  the  requirement  to  develop  good  quality  waveguides  in 
LiNbO^  which  would  confine  both  the  TE  and  TM  modes  equally  well 

Typical  diffused  waveguides  in  LiNbO^  have  a  strong  polari¬ 
zation  dependence.  The  polarization  dependence  of  in-diffused 
waveguides  occurs  because  the  titanium  ions  affect  the  extra¬ 
ordinary  index  more  than  the  ordinary  index.  In  addition, 
unless  careful  procedures  are  introduced,  out-diffusion  causes 
a  refractive  index  profile  for  the  extraordinary  polarization 
separate  from  that  caused  by  in-diffusion.  For  these  reasons 
we  decided  to  investigate  waveguides  fabricated  in  y-cut, 
z-propagating  LiNbO^,  where  both  TE  and  TM  modes  see  an  ordinary 
index . 
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Another  reason  for  being  interested  in  this  orientation 
is  that  it  can  be  used  with  the  electro-optic  effect  to  con¬ 
vert  TE  to  TM  modes,  and  vice  versa.  Either  an  x  or  y-cut 
LiNbO-j  plate  is  used  with  the  channel  waveguide  in-diffused 
along  the  z  direction.  The  electric  field  is  applied  in  the 
plane  of  the  waveguide,  perpendicular  to  the  direction  of 
propagation.  The  input  polarization  can  be  either  TM  or  TE, 
and  the  applied  field  causes  a  rotation  of  the  optic  axes  in 
the  LiNbO^. 

When  this  rotation  is  45°,  a  mode  which  was  originally  pure 

TE  can  be  converted  into  pure  TM,  and  vice  versa.  The  exact 

2 

solution  involves  terms  in  r^  ,  but  these  terms  are  much  less 
than  those  introduced  by  any  alignment  errors  and  can  justi¬ 
fiably  be  ignored.  The  phase  retardation  for  this  case  is 
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The  half-wave  voltage  for  this  orientation , is  only  10-20% 
higher  than  for  the  usual  geometries. 

A  bulk  modulator  of  this  orientation  was  demonstrated  in 
the  experiments  on  optical  bistable  switches  which  were  des¬ 
cribed  in  the  interim  report. 

We  have,  in  principle,  a  design  for  a  multimode  switch, 
and  the  problem  becomes  one  of  making  waveguides  which  guide 
both  TE  and  TM  modes  and  are  thick  enough  to  be  butt-coupled 
to  multimode  fibers.  Out-diffused  waveguides  are  thick  enough. 
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but  act  only  on  the  extra-ordinary  polarization.  Out-diffusion 
increases  the  ordinary  index  but  has  thus  far  not  produced  thick 
enough  waveguides  in  reasonable  times.  One  technique  for  in¬ 
creasing  the  depth  of  the  diffused  waveguides  is  electro¬ 
diffusion,  where,  in  addition  to  the  concentration  gradient 
driving  force,  an  electric  field  is  introduced.  This  method 
has  been  used  to  produce  Cu  in-diffused  waveguides  in  LiTaO^ , 

( Appl .  Phys.  Lett.  2j>  308  (1974))  at  much  lower  temperatures 
and  shorter  times  than  are  usually  required. 

We  have  investigated  titanium  in-diffusion  in  c-propagating 
samples  with  and  without  an  applied  electric  field.  We  have 
results  which  are  inconclusive  and  require  more  experimentation 
before  the  full  understanding  is  developed.  The  work  described 
here  represents  a  progress  report  and  describes  our  status  at 
the  end  of  the  contract. 

Titanium  diffusion  in  the  presence  of  an  electric  field  was 
studied  both  with  electrodes  separate  from  the  samples  and  with 
the  electrodes  in  contact  with  the  samples. 

The  preliminary  measurements  indicated  that  electro-thermal 
diffusion  might  indeed  be  very  satisfactory.  In  an  initial  run, 
waveguides  were  fabricated  by  depositing  500  8  of  Ti  on 
c-propagating  LiNbO^  and  diffusing  in  flowing  C>2  at  1000°C  for 
eight  hours.  The  electro-diffused  sample  was  placed  in  a  field 
of  10  V/mm  during  diffusion,  oriented  in  the  direction  of  the  Ti 
motion.  A  control  sample  was  placed  outside  the  region  of 
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Table  4-1.  Results  of  Titanium  In-diffusion  in  c-propagating 
Waveguides 


Sample 

Number 

Thickness 

Titanium 

8 

Time 

Hours 

Gci  s 

Guided  Mode 
Thickness  (urn) 
TE  TM 

0529  A 

600 

14 

°2 

10 

10 

0713  B 

450 

7 

Ar 

5 

5 

0608  A 

450 

6?5 

Ar 

4 

4 

0608  B 

450 

9 

°2 

5 

5 

0622A 

450 

6^5 

Ar 

10 

10  (near 

cutoff) 

+1 

°2 

— 

— 

+  200 

+  4*5 

Ar 

— 

30  (near 

cutoff) 

0622  B 

450 

85j 

Ar 

— 

20  (near 

cutoff) 

0807  C 

450 

6 

Ar 

— 

20  (near 

cutoff) 

0807  D 

450 

6 

Ar 

— 

20  (near 

cutoff) 

0530  A 

300 

7 

°2 

— 

30  (near 

cutoff) 

0525 

300 

7 

°2 

— 

— 

The 

symbol  —  means 

that  no 

waveguide 

was  observed. 
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the  electric  field.  The  electro-diffused  sample  demonstrated 
good  confinement  of  light  inside  the  waveguide  with  a  high 
N. A.  and  a  guide  thickness  of  70  urn.  The  control  sample 
demonstrated  weak  guiding  in  a  mode  of  30  ym  width  (the 
diffused  region  was  10  ym  thick) ,  with  low  N.A. 

With  these  encouraging  preliminary  results,  experiments 
were  undertaken  to  systematically  study  the  process  of  electro¬ 
diffusion,  and  its  impact  on  the  fabrication  of  waveguides  in 
c-propagating  LiNbO^.  The  results  shown  in  Table  IV-I  were 
obtained. 

The  results  are  somewhat  unreproducible  because  we  had  not 
yet  developed  a  good  method  for  monitoring  the  titanium  film 
thickness,  and  the  numbers  used  here  are  only  approximate.  The 
results  show  that  reasonable  waveguides  can  be  fabricated  with 
film  thickness  of  at  least  450  8.  However,  a  number  of  the 
waveguides  were  right  near  cutoff,  and  the  results  were  somewha4- 
unreproducible .  We  were  surprised  to  find  that  in  many  samples 
there  was  weak  guiding  of  TM  modes,  but  TE  modes  were  not  con¬ 
fined.  This  is  due  to  the  fact  that  the  waveguide  is  anisotropic. 
We  developed  a  theory  to  explain  this  phenomenon,  which  is  in¬ 
cluded  in  the  appendix.  The  result  of  this  study  is  the  conclusion 
that  even  though  the  waveguide  is  c-propagating,  the  TE  and  TM 
modes  do  not  have  the  same  propagation  vectors. 

Further  diffusion  of  waveguides  near  cutoff  caused  all  guiding 
to  disappear.  This  was  also  an  unexpected  result  because. 
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theoretically,  the  longer  the  diffusion  time,  the  farther 
the  waveguide  should  be  from  cutoff.  This  is  shown  theoret¬ 
ically  in  Appendix  II.  Evidently  the  theoretical  models  for 
the  diffused  waveguide  are  not  exact,  when  samples  are  operating 
right  at  cutoff. 

Additional  experiments  were  attempted  with  the  samples 
diffused  in  an  electric  field  and  compared  to  those  diffused 
without  an  electric  field.  The  same  conditions  were  used  as  in 
the  previous  experiments:  450  8  titanium,  6  hour  diffusion  at 
1000°C  in  Argon.  None  of  the  samples  showed  guiding  of  the  TE 
mode.  The  sample  whose  diffusion  occurred  in  the  direction  of 
the  electric  field  showed  no  guiding  in  the  TM  polarization 
either.  A  sample  whose  diffusion  direction  was  counter  to  the 
applied  field,  and  the  samples  which  did  not  have  an  applied 
field  showed  weak  guiding  for  the  TM  polarization,  without  much 
difference  between  them.  For  comparison,  an  X-propagating 
sample  was  studied,  with  a  1000  8  titanium  film  and  2%  hour 
diffusion  in  Argon.  Both  guides  showed  20  pm  TE  modes  and 
10  pm  TM  modes,  with  very  little  difference  between  the  sample 
with  the  applied  field  and  that  without  the  applied  field. 

The  results  of  these  measurements  were  inconclusive  and  did  not 
confirm  the  dramatic  behavior  measured  in  the  first  preliminary 
experiments . 

It  was  decided  to  try  electro-diffusion  with  the  electrode 
contracts  directly  applied  to  the  LiNbO^.  It  was  expected  that 
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the  effects  would  be  bigger  with  this  geometry.  A  sample  was 
prepared  with  titanium  on  the  top  surface,  and  then  both  the 
bottom  and  top  were  overcoated  with  gold.  Gold  wires  were 
attached  in  one  corner  of  the  samples  with  a  thermocompression 
bonder.  The  observations  can  be  summarized  as  follows: 

1.  The  sample  was  heated  to  350°C  for  a  total  of  38  hours  under 
an  electric  field  of  between  10  and  40  V/mm.  It  was  observed 
that  at  temperatures  greater  than  200°C,  the  LiNbO^  sample  began 
to  pass  current.  With  the  voltage  constant  at  20  V/mm,  the  fur¬ 
nace  warmed  up  to  full  temeprature  in  20  minutes,  and  the 
current  rose  from  5  y amp  at  200°C  to  a  maximum  of  32  y  amps,  ten 
minutes  after  the  furnace  reached  its  maximum  temperature.  The 
current  then  began  to  decrease  to  about  5  y  amps  about  8  0  minutes 
after  the  oven  reached  its  maximum  value.  After  a  long  period 
of  steady  state,  the  electric  field  was  suddenly  switched  to 

40  V/mm.  The  current  rose  slowly  again  to  about  30 y amps, 
reaching  a  peak  40  minutes  after  the  change  of  electric  field. 
After  saturating  at  its  maximum  value,  the  current  slowly  de¬ 
creased  to  17  y  amps,  170  minutes  after  the  voltage  change. 

At  the  temperature  of  350°C,  when  the  electric  field  was 
suddenly  switched  from  0  to  a  final  value  of  20  V/mm,  the  cur¬ 
rent  responded  immediately  without  any  delay,  to  a  new  peak 
value. 

2.  During  the  experiment.  Ohm’s  law  did  not  hold;  i.e.  the 
resistivity  appeared  to  be  voltage  dependent.  These  results 
are  shown  in  Table  4-II.  The  cause  is  not  understood. 

3.  After  electro-diffusion  at  350°C,  there  was  a  crack  across 
the  substrate.  In  addition,  there  was  no  evidence  of  a  waveguide. 


Table  4-II.  Resistivity  Measured  at  300°C  and  600°C 

Resistivity 


Voltage 

Calculated  from 

Measured 

Temp. 

Across  the  Sample  Current 

Ohm ' s  Law 

Resistivity 

°C 

V 

mA 

cm 

cm 

300 

12 

2.5 

X 

10-3 

6.0  x 

io7 

18 

4.3 

X 

10~3 

5.2  x 

io7 

40 

18.1 

X 

io“3 

2.8  x 

107 

600 

6 

14 

5,4  x 

io3 

2.5  x  104 
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The  sample  came  out  black  in  the  vicinity  of  the  wire  bond. 
Studies  of  this  darkened  area  on  the  Xray  dispersive  SEM  in¬ 
dicated  Si  and  Mg  contamination.  It  is  not  clear  how  these 
metals  got  in  from  the  bond,  but  it  was  decided  to  redesign 
the  electrode  arrangement.  No  waveguiding  was  observed  with 
this  sample;  it  was  too  dark. 

A  new  electrode  arrangement  was  developed,  whereby  the 
sample  was  mounted  between  two  platinum  plates.  With  this 
arrangement  no  darkening  occurred.  However,  cracks  still 
appeared.  Only  one  sample  cut  of  four  was  not  cracked.  That 
sample  also  was  the  only  sample  that  did  not  pass  current 
during  the  diffusion  process.  This  leaves  the  experiments  at 
a  very  interesting  stage.  Perhaps  the  current  is  caused  by  the 
crack  and  when  there  is  no  cracking,  there  is  no  current. 
Perhaps,  instead,  there  was  not  good  electrical  contact  for 
that  one  sample.  That  sample  was  the  only  one  that  could  be 
tested  for  waveguiding.  None  was  observed.  However,  until  we 
demonstrate  clearly  that  there  was  electrical  contact  with  the 
sample  during  diffusion,  we  cannot  conclude  that  electro-thermal 
diffusion  was  not  successful. 

We  suggest  that  it  would  be  useful  to  continue  experiments 
on  electro-thermal  diffusion  to  develop  the  conditions  under 
which  it  works.  Since  the  effect  has  been  reported  in  LiTaO^, 
we  expect  that  we  should  be  able  to  develop  it  for  LiNbO^. 

One  advantage  of  electro-thermal  diffusion  is  the  potential 
for  driving  the  waveguide  below  the  LiNbO-j  surface.  This  has 
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two  distinct  advantages:  first,  it  reduces  losses  due  to 
surface  scatter;  second,  it  causes  a  symmetric  refractive  index 
profile  and  results  in  a  mode  profile  that  more  closely  matches 
that  of  the  optical  fiber.  The  theoretical  analysis  of  the 
effect  of  the  applied  field  on  the  mode  profile  is  shown  in 
Appendix  III. 

In  conclusion,  the  experiments  on  electro-thermal  diffusion 
in  LiNbO^  do  not  yet  show  any  dramatic  evidence  for  multimode 
guides  which  confine  both  TE  and  TM  modes.  More  experiments 
should  be  done  to  determine  whether  or  not  this  technique  will 
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APPENDIX  4-1  TM-Selective  Guiding  in  c-propagating  LiNbO^ 

TM-selective  guiding  is  observed  near  cut-off.  This  Appendix 
contains  a  theoretical  explanation  using  the  model  of  a  step-index 
planar  waveguide  in  an  anisotropic  medium,  with  the  same  symmetry 
as  LiNbO-j. 

For  simplicity  we  choose  a  weakly  guiding  model,  assuming 
that  the  refractive  index  components  of  the  guiding  layer  are  only 
slightly  larger  than  those  of  the  bulk. 

Abbreviations  to  be  used  are  as  follows: 

nQ  =  ordinary  refractive  index  of  bulk 

ne  =  extraordinary  refractive  index  of  bulk 

NQ  =  ordinary  refractive  index  of  waveguide 

Ne  =  extraordinary  refractive  index  of  waveguide 

k  =  wave  vector 

6  =  displacement  vector 

8  =  angle  between  the  optic  axis  (propagation  direction) 

and  the  direction  of  the  ray  propagating  in  the  waveguide. 

I 

8  =  angle  between  the  optic  axis  (propagation  direction)  and 

the  direction  of  the  ray  propagating  in  the  bulk. 

B  =  propagation  constant  of  a  guide  mode 

N  =  refractive  index  of  an  extraordinary  ray  in  the  guide 

n  =  refractive  index  of  an  extraordinary  ray  in  the  bulk, 

b  =  waveguide  thickness 

kQ  =  wave  number  in  free  space 


Because  of  the  anisotropy,  it  is  required  to  talk  about  index 
ellipsoids  rather  than  refractive  indices.  In  this  theory  we 
employ  two  refractive  index  ellipsoids,  in  the  bulk  and  waveguide 
respectively.  The  waveguide  ellipsoid  is  slightly  enlarged  com¬ 
pared  with  that  of  the  bulk  due  to  the  Ti-indif fusion  and  Li 
out-diffusion.  In  Fig. 41-1  the  two  ellipsoids  are  shown.  The  index 
increase  is  shown  to  be  much  larger  along  the  c  axis  than  perpen¬ 
dicular  to  it.  This  note  will  describe  why  this  is  so. 

In  the  ray  optics  picture  of  propagation  in  a  waveguide,  the 
guided  modes  can  be  considered  to  be  bouncing  back  and  forth  between 
the  waveguide  walls  at  an  angle  9.  Those  rays  which  are  polarized 
TE,  are  pure  ordinary  rays  and  do  not  experience  the  anisotropy 
of  the  waveguide.  The  TM  modes,  however,  have  a  polarization  which 
lies  in  the  plane  of  incidence  of  the  rays,  and  have  a  component 
of  polarization  which  lies  along  the  c-axis.  This  is  shown  in 
Fig. 41-2. 

We  know  from  birefringence  theory  that  the  TM-polarized  ray 
shown  in  Fig.  2,  traveling  at  an  angle  9,  will  experience  a 
refractive  index  given  by 

N  N 

N(9)  =  - * - * - - — - 5 - j-  (1-1) 

(Ng  cos  0  +  Nq  sin  0)^ 

The  light  traveling  in  the  substrate,  which  has  an  exponential 
profile,  can  be  considered  to  be  composed  of  rays  traveling  at 
an  imaginary  angle,  0'.  This  makes  it  possible  to  write 
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(1-2) 


Ncos©  =  n  cos  ©'  =  6/k 


with 


n  (©' )  = 


n  n 
e  o 

(n  cosz©'  +  n  sin  ^©') 
e  o 


(1-3) 


N  and  n  are  written  in  terms  of  3 


2  1  M  2M  2  2  ..  2  3 
N  =  - ~  N  N  _  -  (N  -  N  - » 

N  2  e  °  e  °  k  - 

o  o 


(1-4) 


oi  2  2  ,  2  2,  . 

2  1  n  n  -  (n  -  n  )  3 

n  =  — =-  e  o  e  o  — 

n  k 

o  o 


(1-5) 


We  can  write  the  waveguide  equation  in  the  bulk 


-C  +  (n2k  2  -  B2)  PI  =  0 

UJ 


(1-6) 


and  in  the  waveguide 


9 2  +  (N2k  2  -  B2) 

2  ° 

9y 


IB 


(1-7) 
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The  condition  for  guided  waves  is 


2  2  2  2  2  2 
(n  Kq  -  8  )  <0  and  (!Tk0  -  g  )  >0 

If  we  use  abbreviations 

v  =  (8  -  n  )^b  and  u  =  (*rkQ  -  8  Pb 


we  write 


2  2 
u  >  and  u  >  0 . 


For  TE  inodes  the  relationship  between  u  and  v  is  the  same  as 
for  isotropic  guides: 


where 


2  2  „2 

u  +  v  =  V 


V  -  «02  -  "o2>  V 


(1-8) 


For  TM  modes,  however,  u  and  v  satisfy  the  equation  for  an 
ellipse,  which  does  not  appear  in  the  isotropic  guide: 


u  +  v  -  1 

+  _  i 

A  B 


(1-9) 


where 


A  =  ^  (N  2  -  n  2)  **  k  b 

N  O  O  O 

O 


and 


B  =  ~ •  (N  2  -  n  2) **  k  b 
n_  o  o  o 
o 
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To  obtain  the  eigenvalue  equation,  we  apply  the  same  analysis 
as  applied  for  isotropic  waveguides,  since  all  the  procedure  in  the 
boundary  condition  matching  between  the  waveguide  and  the  bulk  must 
be  followed  exactly  in  the  same  fashion.  Thus  we  can  use  the  result 
of  a  weakly  guiding  isotropic  step  guide  to  get  for  both  TE  and  TM 
modes 

v  =-u  cot  u  (1-10) 


We  have  confirmed  that  the  quasi-ray  analysis  employed  coincides 
with  the  results  of  the  field  treatment  done  by  Yamamoto  et.al.^ 
The  eigenvalues  for  B  will  be  determined  by  the  crossing 
points  of  Eqs.  (9)  and  (10) for  TM  mode  and  Eqs.  (8)  and  (10)  for  TE 
mode.  In  Fig.  3  the  characteristic  Eq.  (10)  and  Eqs.  (8)  and  (9) 
are  plotted.  From  this  we  can  readily  determine  the  three  different 
cases  of  guiding: 


With  the  extraordinary  index  increase  Ane  and  the  ordinary  index 
increase  AnQ  in  the  guiding  layer  we  can  relate  the  refractive 
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indices  of  the  bulk  and  the  waveguide  by 

Ne  =  n  +  An  and  N  =  n  +  An 
e  e  o  o  o 

For  LiNbO^  (n  =  2.29  and  ng  =  2.20  at  0.63  ym)  the  condition  of 

TM  mode  -  selective  guiding  is  determined  to  be 

(An  -  An  )  >0.09. 
e  o 

We  should  point  out  that  the  step  guide  analysis  is  not  good 

for  the  quantitative  analysis  of  a  real  graded  index  guide  and 

therefore  we  cannot  judge  this  theory  from  the  value  calculated 

(2) 

with  this  model.  However,  an  estimate  of  An  -  An  =0.03 

e  o 

is  typical  for  Ti-indif fused  guides.  Thus  we  consider  that  the 
above  theory  is  undoubtedly  the  explanation  for  the  TM-only 
guides  which  we  observed. 
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APPENDIX  4-II.  Waveguide  Cutoff  Condition 
Dependence  on  Diffusion  Time 

Let  us  consider  diffusion  from  a  finite  source  of  diffusant, 
such  as  a  thin  film  of  LiNbO^*  The  refractive  index  will  have  its 
peak  value  at  the  surface,  where  it  will  be  n  larger  than  its 

value  in  the  bulk,  N^.  As  the  diffusion  time  gets  longer,  the 

% 

surface  value  of  An  decreases  by  t  ,  while  the  diffusion  profile 

r 

spreads,  increasing  the  diffusion  depth  by  t  .  Thus,  while  the 
index  step  decreases,  the  effective  waveguide  thickness  increases. 
This  contradictory  behavior  of  the  two  waveguide  parameters  does 
not  allow  us  to  make  a  quick  judgment  as  to  whether  increasing  the 
diffusion  time  will  increase  or  decrease  the  waveguide  cutoff.  In 
this  appendix  we  examine  the  cut-off  condition  for  a  diffused  guide 
index  profile. 

It  is  well  known  that  for  a  sufficiently  long  time  diffusion  of 
a  limited  amount  of  diffusion  source,  the  diffused  concentration 
profile  follows  the  Gaussian  distribution.  We  can  reasonably 
assume  the  refractive  index  change  will  follow  the  same  profile. 

In  order  to  judge  whether  a  sample  forms  a  guiding  layer,  it  is 
necessary  only  to  determine  the  cut-off  condition  for  a  single  mode. 
The  lowest  order  mode  cut-off  condition  is  ^ 


An 


dx  -  3  tt/4 


(1) 


where  k  is  the  free  space  wave  number,  x  is  the  position  variable 

measured  from  the  surface  and  b  is  the  —  width  of  the  Gaussian 

e 

index  profile.  Equation  (1)  can  be  calculated  to  give 


kQ  b  2  n^  An  =  3  tt/4 


(2) 


For  finite  source  diffusion,  the  quantities  n  and  b  are 
related  by 


An  _2  dn 

t/b  tt  dc 


(3) 


where  t  is  the  deposited  metal  film  thickness, a  is  the  metal  atom 

concentration  of  the  film  and  dn/dc  is  the  rate  of  change  of  the 

refractive  index  with  diffusant  concentration.  The  quantitity 

(2) 

dn/dc  has  been  shown  to  be  constant  by  Schmidt  and  Kaminow  as 

well  as  Naitoh  et.al.  ^  Two  sets  of  measurements  give  almost 

-23  3 

the  same  value,  dne/dc  =  1.6  x  10  cm  ,  for  the  extraordinary 
index  change  due  to  Ti-  diffusion.  Now,  if  we  assume  dn/dc  to 
be  constant,  substitution  of  Eq.  (3)  into  Eq.  (2)  gives 

kQ  2  n^  t  (dn/dc) ab  =  j  n  .  (4) 
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In  other  words,  the  film  thickness  required  to  ensure  confinement 
of  at  least  one  mode  is 


(5) 


k 

Since  b  is  proportional  to  t  ,  a  longer  diffusion  time  decreases 
the  film  thickness  required  to  produce  a  waveguide. 

For  a  given  deposition  thickness  of  diffusant,  the  longer  the 
diffusion  time,  the  better  confined  the  mode  should  be.  This 
theoretical  result  does  not  agree  with  our  experiments.  This  is  not 
totally  unexpected,  however,  because  the  left  hand  side  of  Eq.  (4) 
is  a  slowly  varying  function  of  time.  A  little  discrepancy  between 
the  assumptions  made  in  this  theory  and  the  real  diffusion  para¬ 
meters  could  revert  the  increasing  time  dependence.  After  a  long 
diffusion  time,  the  refractive  index  profile  may  deviate  substan¬ 
tially  from  the  theoretical  Gaussian  profile.  In  addition  to  Ti 
in-diffusion,  there  are  three  other  diffusion  processes  which  may 
be  participating  in  forming  the  refractive  index  profile.  Those 
are  suggested  as  Li20  out-diffusion,  Nb02  diffusion  and  oxygen 
revolution  with  Nb^+  Nb^+  reduction/^  Near  1000°C,  Li20  and 
oxygen  diffusions  are  significant  but  NbQ2  diffusion  is  not, 
because  of  its  extremely  low  vapor  pressure.  Li20  out-diffusion 
is  well  known  to  affect  only  the  extraordinary  index.  It  is  not 
known  to  us  how  much  the  oxygen  diffusion  affects  the  refractive 
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index  profile.  However,  we  expect  that  it  could  have  a  large 

effect  because  an  oxygen  ion  has  bigger  electronic  polarizability 

4+  ( 4 ) 

3.88  in  c.g.s.  unit  than  0.19  of  Ti  .  Thus  the  resultant 

refractive  index  profile  after  diffusion  in  LiNbO^  will  be  deter¬ 
mined  as  well  as  the  superposition  of  diffused  Ti  contribution 
and  the  side  effects  due  to  other  three  contributions  and  may 
deviate  from  the  Gaussian  profile  based  on  a  single  ideal  diffusion 
theory. 


REFERENCES  I-II 


1.  H.  Naitoh,  M.  Numoshita,  and  T.  Nakayama,  "Mode  Control  of 
Ti-diffused  LiNbO^  Slab  Optical  Waveguide",  App.  Opt.  16 , 

2546  (1977). 

2.  R.V.  Schmidt  and  I.P.  Kaminow,  "Metal-diffused  Optical 
Waveguides  in  LiNbO^",  App.  Phys.  Lett.  2jj,  458  (1974). 

3.  I.P.  Kaminow  and  J.tR.  Carruthers,  "Optical  Waveguiding  Layers 
in  LiNb03  and  LiTa03",  Appl.  Phys.  Lett.  22,  326  (1973). 

4.  C.  Kitlel,  "Introduction  to  Solid  State  Physics",  5th  ed.  , 

John  Wiley  &  Sons,  Inc. ,  New  York/London/Sydney/Toronto, 

411  (1976). 


APPENDIX  4-III 


Previous  experiments  on  buried  waveguides  have  used  a  double 
diffusion  technique,  following  an  in-diffusion  of  one  species  to 
raise  the  index,  with  another  to  lower  the  index  closer  to  the 
surface.  Theoretically  a  single  diffusion  with  the  help  of  an 
electric  field  can  make  a  buried  wave  guide.  The  diffusion 
equation  under  an  external  field  is 


3c 
3 1 


(III-l) 


where  C  is  the  diffusant  concentration  in  the  substrate,  x  is  the 
coordinate  measured  from  the  surface  of  the  sample,  D  is  the 
diffusion  constant,  y  is  the  mobility  of  the  diffusant,  and  E  is 
the  applied  external  field. 

It  is  well  known  that  the  diffusion  profile  without  an 
external  field  becomes  Gaussian  for  a  long  diffusion  time  from  a 
limited  source.  The  solution  of  the  diffusion  equation  without 
a  field  is: 

2C  9  , 

c.  =  —  exp  (  -  xVb  )  (II I— 2 ) 

firb 


where  b  is  the  diffusion  depth  and  cQ  is  the  number  of  metal  ions 
per  unit  area  which  are  available  for  the  diffusion  (CQ  metal- 
thickness)  . 

In  order  to  solve  equation  III-l,  take  a  coordinate  trans¬ 
form  such  that  x  x'  +  yt' 


The  diffusion  equation  then  has  a  very  simple  form: 


3c  (  x' ,  t' )  _  _  3  C  (  x' , t) 
9t  ‘  3x2 


Immediately  we  get  a  solution  of  the  form 

2c 

C  (x't)  =  A  exp  (-x12/b2) 

tfTb 

2C  ,  ■> 

C  (  x,t)  =  A  — —  exp  {  -  (  x  -  E]it )  /b  ) 
mb 


The  new  concentration  represents  a  buried  profile  below  the 
surface  by  a  distance  XQ  =  Eut.  The  constant.  A,  in  the  above 
equation  is  determined  by  requiring  conservation  of  the  amount 
of  diffusant.  Clearly  the  waveguide  can  be  buried  most  effec¬ 
tively  if  the  diffusant  has  a  relatively  high  electrical 
mobility  in  the  substrate. 

This  guided  mode  profile  matches  very  well  the  symmetric 
profiles  from  optical  fibers.  Electrical  burying  of  the  wave¬ 
guide  below  the  surface  would  lower  the  insertion  loss  of 
LiNb03  switches  used  with  fiber  optics. 
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